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STRESSES  AND  DISPLACEMENTS 

IN  TWO,  THREE  AND  FOUR  LAYERED  SYSTEMS 

WITH  FIXED  BOTTOM 


Introduction 

This  report  deals  with  the  mathematical  aspects  required  for  the  establishment 
of  a  computer  program  able  to  calculate  all  stresses  and  displacements  in  two, 
three  and  four  layered  systems. 

The  materials  of  the  different  layers  may  be  isotropic  or  cross-anisotropic. 

The  interface  conditions  cover  all  the  cases  from  full  friction  to  full  slip  included 
partial  friction. 

The  bottom  or  the  last  layer  is  considered  to  be  fixed  (no  vertical  deflections). 

The  loads  can  either  be  flexible  either  rigid. 

The  report  is  based  on  : 

-  existing  material  :  isotropic  multilayer  theory  (BURMISTER,  1943) 
and  anisotropic  multilayer  theory  (VAN  CAUWELAERT,  1983). 

-  original  research  work  :  interface  conditions  (fixed  bottom,  partial  friction)  and 
satisfactory  convergency,  thus  complete  accuracy,  at  the  surface  and  in  the  first 
layer  of  the  system,  rigid  load  boundary  condition. 

This  report  contains  three  parts  : 

Part  1  :  a  theoretical  outline  wherin  the  basic  equations  are  given,  the  specific 
boundary  conditions  discussed  and  the  particular  numeral  problems  related 
to  the  accuracy  at  the  surface  and  in  the  first  layer,  solved. 

Part  2  :  the  general  mathematical  analysis  of  the  chosen  numerical  solution  and 
a  description  of  the  programs  and  their  utilization. 

Part  3  :  (appendices)  the  detained  mathematical  and  algebrical  analysiseas  for  the 
different  considered  cases  :  isotropic  or  anisotropic,  full  slip  or  not,  flexible 
or  rigid  load. 

The  programs  are  written  in  FORTRAN  77  and  run  on  IBM  PC  or  all  other  compa¬ 
tible  equipment. 


m. 


JL 
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PART  1.  THE  MULTILAYER  SOLUTION 

1.1  THE  BASIC  EQUATIONS  AND  FUNDAMENTAL  HYPOTHESISES 

The  itrtiau  and  displacements  In  a  multilayered  system  of  homogeneous  and 
isotropic  layers  subjected  to  a  uniform  vertical  load  applied  over  a  circular 
area  are  obtained  from  the  following  stress  function: 

J  (mr).F(m) 

♦  -  fo  ~ — ; - IV“  -  Bi«““  ♦  lCt«“  -  ‘V-**1  d“ 

Resulting  stresses  and  displacements  are  given  by: 

,•  2  am  2  it 

°z  *  /©  J0(«)-^(*)lAim  e  +  B^m  e 

-  C^m(l  -  2nt  -  m*)e**  ♦  D^md  -  2^  ♦  m)*"** ]  dm 

ac  m  ~  fQ  J0(mr).F(m)  lA^e"*  ♦  B^c”” 

+  Cjmd  ♦  2ut  ♦  m*)e"*  -  Dtm(l  ♦  2ui  -  maje”**]  dm 


J^mrj.FCm) 

— Si - l\"  *“ 


♦  C^mC 1  ♦  m»)e"*  -  D^mC 1  -  ms)e~** 1  dm 

<?q  •  ~  J0  J0(mr).F(m)lCim.c"*  -  D^m.e”** J .2^ .dm 
^J. (mr).F(m)  , 

-  /0 - — - lA^m  e  ♦  B^m  e"**  ♦  Ctm(l  ♦  maje"* 

-  0^(1  -  rnsje”**  j  dm 
*  ”  /p  Jj(«r).F(m)  lA^m^e"*  -  B^e"** 


i  •  radius  of  cha  uniformly  load ad  circular  area 
p  -  intensity  of  tha  uniform  vertical  load 

r  •  horizontal  distance  from  tha  axis  in  a  cylindrical  coordinate 

system 

z  ■  depth 

0^  “  vertical  stress 

o  “  horizontal  radial  stress 
t 

o  -  circumferential  stress 

-  shear  stress 

v  ■  vertical  deflection 
u  •  radial  (horizontal)  displacement 

-  Young  modulus  of  a  given  layer 

•  Poisson's  ratio  of  a  given  layer 

•  unknown  parameters  to  be  determined  by  the  boundary  conditions 

J  •  Bessel  function  of  the  first  kind  of  order  zero 
o 


a  ■  integration  parameter 
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la  Che  case  of  a  cross-anisotropic  body  Che  screes  function  is: 


mo  -ms 

t  *  /"  Jo(ar).F(«)/m  [A^e®*  -  B^"®*  +  (^e  *2  -  D^e  i2]  dm 

Thf*  stress  function  differs  fundamentally  from  the  isotropic  case  such  that 
the  two  casea  must  be  handled  different. 


The  stresses  and  displacements  are  given  by 

°z  “  Jo  J0<Br)*F<*)  +  ut)  <A1m2e"*  +  Bim2c-nZ) 


+  n^Cn^  +  w^)  (Cis^n^e"’*i,,*  +  D^S^m^e  im2)]  dm 
“  JQ  JQ(mr).F(m)  (n^l  ♦  ui}  (AtmV“  +  B^n^e"®*) 


n.  (n.  -  p .  )  _ 

+  Tt~t  (cisi"  eV“  +  Di\®  d“ 


.  J.(»r)F(m) 

[0  - — -  Qt(l  ♦  w1)[Aim  e  +  e 


2  *l“  2  -*i“ 

+  C^s^m  e  +  D^s^m  e  1  dm 

•  “•  ms  n  p  (  l-n  ) 

-  C  I'iV  *»lV  I  «/-  ,  * 


J  (mr).F(m)  , 

*  fa  _ . 


}0  ST 


2  ms  2  -ms  2  *i“ 

jAtm  e  ♦  B^m  e  ♦  C^m  e 


♦  1  J  *»1(l+ut)dm 


tfl  m  ~  I*  Jj<«r).F(m)  lOjU  ♦  (A^e**  -  »1m2e*®*) 


and  that  v,  Poisson's  ratio  in  ths  horisontal  plane,  is  related  to  Poisson's 
ratio  in  the  vertical  plana  by  (EFTIMIE,  1973;  VAN  CAUWELAERT,  1983): 
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vi  ’  “i/ni 


1.2.  THE  BOUNDARY  CONDITIONS 

1.2.1  The  Surface  Conditions 

The  surface  conditions  are  expressed  by  an  adequate  value  of  the  Kernel 
F(m)  in  the  stress  function. 

In  the  case  of  a  uniform  distributed  load  over  a  flexible  area  the  kernel 


F(m)  *  paJ^(ma) 

where, 

I 

Jj  *  Bessel  function  of  the  first  kind  of  order  one. 

Indeed  the  expression  for  the  vertical  stress  is  then  at  the  surface. 

-00 

\  o  *  paf  J  (mr)  J.(ma)  dm 

z  '  o  o  1 

•  p  for  r  <  a 

■  p/2  for  r  *  a 

■  o  for  r  >  a 

In  the  case  of  a  load  distributed  over  a  rigid  area  the  kernel  is 

F(m)  -  ^  sin  (am) 


JL... 


.a . . 
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at  the  surface 

cr  “  ¥r~  f  J  (mr)  sin  (am)  dm 
z  zoo 

£a  ,  2  2-1/2 

-  y-  (a  -r  )  for  r  <  a 

■  0  for  r  >  a 

The  total  load 


P  =  fa  o  rdrdQ 
‘  o  1  o  z 


2 

irpa 


is  equal  to  the  applied  load. 

Also  at  the  surface,  the  deflection  is 
2 

(1-U.)  J  (m)  sin  (am) 

w  -  — pa r  - -  dm 

E .  J  o  m 


2U-v h 

- - p - for  r  <  a 

E1  * 

The  deflection  is  constant  under  the  load. 


1.2.2  The  Interface  Conditions. 


Let  as  consider  an  n-layered  system,  consisting  of  (n-1)  layers  each  hav¬ 
ing  a  finite  thickness  supported  by  a  semi-infinite  body.  There  is  associated 
with  each  layer  a  stress  function  (A^  D^)  with  4  unknown  parameters 

such  that  the  total  of  unknown  parameters  is  4n.  Two  of  the  parameters  depend 
on  the  loading  surface  conditions. 


o 

z 


f(p)  for  r<  a 
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At  infinite  depth  the  stresses  and  displacements  must  vanish  and  thus  A  and 

n 

Cq  ■  o.  As  a  result,  there  are  4  (n-1)  parameters  to  be  determined  with  4 
conditions  at  each  interface.  The  four  interface  conditions  are  addressed  by 
imposing  the  conditions  that  the  layers  maintain  contact  and  that  the  vertical 
stresses  (o z  ),  shear  stresses  (x  )  and  vertical  displacements  (w)  at  the 
bottom  of  each  layer  and  at  the  top  of  the  underlying  layer  are  equal.  The 
fourth  interface  condition,  horizontal  displacements  (u),  depends  on  the  rela¬ 
tive  adhesion  at  the  interface  between  the  layers.  The  two  extremes  of  adhen- 
sion  are 


a  full  continuity,  expressed  by  setting  the  horizontal  displacement  (u)  on 
each  side  of  the  interface  equal. 

®  frictionless  interface  expressed  by  considering  the  interface  as  a  prin¬ 
cipal  plan  which  results  in  the  shear  stresses  equaling  zero. 

Partial  adhesion  can  be  expressed  by: 


with 


K  e  [O.inf] 

When  K  ■  1,  one  has  full  continuity 
K  *  1,  one  has  partial  continity 

Zero  friction  then  becomes  a  separate  case,  for  which  another  program  has  to 
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be  written. 

1.2.3  The  fixed  bottom  condition 

The  boundary  conditions  discussed  in  the  previous  paragraph  indicate  that  the 
last  layer  of  the  multilayer  is  considered  as  a  semi-infinite  body.  One  can 
also  consider  the  case  of  a  multilayer  system  resting  on  a  rigid  body,  such 
that  vertical  displacements  vanish  at  the  contact  face  with  the  rigid  body. 
This  problem  shall  be  designated  as  the  fixed  bottom  condition. 

In  addition,  through  the  general  solutions  of  the  compatibility  equation 

in  multilayer  layer  theory,  we  can  get  w=o  at  the  desired  depth  and  determine 

the  corresponding  values  of  the  parameters  A  ,B  ,C  ,D  .  In  the  isotropic 

n  n  n  n 

case,  the  parameter  C  must  then  be  put  equal  to  zero,  while  the  parameter  A 

n  11 

is  determined  by  condition  w=“0.  In  the  anisotropic  case,  either  the  parameter 

A  must  be  zero  when  s<l,  or  the  parameter  C  must  be  zero  when  s>l  Using  cri- 
n  n 

teria  that  a  minimum  influence  on  the  surface  deflection  is  desireable  and  the 

fact  that  conditions  u  *  o  and  T  =  o  have  less  physical  sense,  we  retain  the 

rz 

condition  w  ■  o  as  the  most  reasonable  rigid  bottom  condition.  An  added  bene¬ 
fit  of  this  selection  is  that  the  condition  w=o  is  the  easiest  to  account  for 
mathematically . 

1.2.4  Influence  of  a^  fixed  bottom  condition  on  the  numerical  computation. 

The  deflection  at  the  surface  of  a  two  layered  system,  without  fixed  bottom, 
was  solved  by  BURMISTER  (1945) 

2(1 -Uj)  ^  J0(mr> .Jj(ma) 
w  "  "Pa  — -  /  o  - m - 
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-2mh  — 4mh 

1  +  4Kmhe  -  KLe 


1  -  (L  +  K  + 


4km2h2)« 


-2mh 


+  KLe 


-4mh 


dm 


where 


K  = 


1  -  n 


(3-4p2)  -  n(3-4Wl) 


1  +  n(3-4u1) 


(3-4u2)  +  n 


E2  +  «*!> 
E1  ^  +  y2) 


which,  to  avoid  convergency  problems  during  the  numerical  computation  (see  art 
1.3.3),  is  written  as: 


2(1  -  u)  J  (mr).J  (ma) 

1  (»  o  1 

w  *  “pa - v -  L  - - - dm 


‘1 


m 


-pa 


.  2)  J  (mr).J.(ma) 

2(  1  -  U  r®  O _ 1 _ 


El  jO 


m 


2  2  -2mh  -4mh| 

(LUt  4Kmh  +  4 Km  h  )e  -  2 KLe  ■ 


1  -  (L  +  K  +  4Km2h2)e”2mh  +  KLe_4mh 


dm 


The  deflection  with  depth  on  the  axis  of  the  loaded  area  is  computed  (r  -  0, 

J  (mr)  *  1)  with  E,  -  1,000,  E  ■  100, p  ■  p0  -  0.5  ,  a  -  10,  h  -  10,  p  *  1. 

o  1  Z  I  z 

The  value  of  is  obtained  analytically  (WATSON,  1966) 
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2(1  -  UJ)  „  J^ma)  2(1  -  uh 

w,  »  -pa - e -  JQ  — j -  dm  -  -pa  - g -  =  -0.015 


'1 


The  value  of  is  obtained  by  numerical  integration  for  different  values  of 
the  integration  parameter  m: 


m 

w 

2 

w 

5 

t 

f 

0.5 

-  0.0582 

-  0.0732 

» 

> 

0.1 

-  0.0580 

-  0.0730 

0.05 

-  0.0580 

-  0.0730 

0.02 

-  0.0580 

-  0.0730 

« 

0.01 

-  0.0580 

-  0.0730 

From  this  data,  the  deflection  is  correct  for  values  of 

For  the 

case  of  a  rigid  bottom 

at  a  depth  H 

the  value  of  the 

surface 

is  given,  with  u  *  0.5, 

by 

_  1.5pa 

J^ma) 

,  -2mH 

1  -  e 

E  ■'o 

m  1  +  (1  +  2mH)e_2mh 

In  this 

expression,  for  m  *  0, 

the  value  of 

the  integrand  is 

dm 


while  in  the  previous  case  the  value  was : 


-pa 


2(1  -  up  , 


L  +  K  -  2KL 


*  2  *  1  -  L  -  K  -  KL 


i 

I 
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By  splitting  the  integral  the  expression  becomes: 


J  (ma) 
1.5par®  1 

w  * - s1— I  - dm 

E  *  n  in 


+  UHf 

E  Jo 


J  j (ma ) 


2(  1  +  mH)e  _ 
1  +  (1  +  2mH).e' 


With  p  -  1  ,  a  ■  10  and  E  -  1,000,  the  value  of  ■  -  0.015;  the  values  of  w^ 
as  a  function  of  m,  are  computed 


-  with  H  -  3a 


m 

w2 

w 

0.5 

+  0.0036 

-  0.0114 

0.1 

+  0.0033 

-  0.0117 

0.05 

+  0.0032 

-  0.0118 

0.02 

+  0.0032 

-  0.0118 

-  with  H  -  1,000a 


0.5 

+ 

0.0833 

+  0.0683 

0.1 

+ 

0.0167 

+  0.0017 

0.05 

+ 

0.0083 

-  0.0067 

0.02 

+ 

0.0033 

-  0.0117 

0.01 

+ 

0.0017 

-  0.0133 

For  the  case  of  a  rigid  bottom  at  H  -  3a,  the  deflection  is  correct  for  values 
of  m  ■  0.05  or  smaller  but  for  the  case  of  a  rigid  bottom  H  *  1,000a,  the 
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deflection  ie  not  correct  for  all  the  values  a.  In  the  latter  case  (H  - 
l 000a),  the  problea  is  analogious  to  the  semi-infinite  layer  case  where  w  - 
-0.01. 

For  a  *  0,  the  function  f(a),  in  1.5  pa/E  /  F(a)  da,  is  equal  ?ero:  but  at  a 

o 

value  of  a  •  CH-,  the  function  is  equal  to  its  seai-inf inite  body  value. 

For  coaparison,  the  function  is  given  below  for  various  values  of  a  for  a 
seal-infinite  body  and  for  a  deep  rigid  bottoa. 


a 

f(a)  seal-infinite 

f(a)  fixed  bottoa 

0.00 

0.5 

0 

0.01 

0.499999 

0.499999 

0.02 

0.499975 

0.499975 

0.05 

0.499844 

0.499844 

0.10 

0.499375 

0.499375 

In  choosing  a  value  too  high  for  a,  an  important  part  of  the  deflection  is 
neglected.  The  error  is  significant  for  the  deep  rigid  bottoa  condition.  As 
a  result,  the  rigid  bottoa  case  aust  be  handled  carefully,  while  an  average 
value  of  a  ■  0.1  gives  correct  results  for  the  seai-inf inite  case. 

1.3.  PARTICULAR  NUMERICAL  PROBLEMS 

Several  numerical  probleas  arise  when  accuracy  is  desired,  whatever  the 


location  at  which  stresses  and  displaceaents  are  to  be  computed.  In  all  cases 
the  probleas  of  accuracy  can  only  be  solved  where  the  relations  for  stresses 
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and  displacements  ara  aval  labia,  although  partially,  in  doaad  to  cm.  Caraful 
at t ant loo  to  a pacific  terns  within  tha  gaoaral  expressions  can  halp  to  improve 
tha  accuracy  of  calculationa. 

1.3.1.  Tha  full  alip  lntarfaca  condition 

As  was  pointed  out  tha  xaro  friction  or  frictional  ess  interface  condition 
ia  a  separata  case  that  will  be  addressed.  Tha  value  of  any  atraaa  or  dis¬ 
placement  is  obtained  froa  one  of  tha  above  aentioned  relations.  Let  us  con¬ 
sider,  for  axaapla,  tha  vertical  stress  in  tha  i-th  layer  of  an  isotropic 
layer: 

■  P*  /c  J0(«r).J  j(aa)  (A^e**  ♦  Bta2a”“  -  Cta<  l-2wt-aa)e“ 

♦  U^at  l-Zuj^+asJa-"*  1  da 

Solution  of  this  relation  by  nuaarical  integration  requires  a  value  of  a  froa 
o  to  a  value  high  enough  to  ensure  convergency.  In  the  case  of  n  layers  and 
rigid  bottoa  the  paraswters  A^.B^.C^  and  aust  be  deterained  for  each  value 
of  a  baaed  on  boundary  conditions  and  a  systea  of  (An  -  l)  equations  with  (An 

-  1)  unknowns. 

Early  solutions  of  this  problea  involved  inverting  the  awtrix  of  the  (An 

-  1)  unknowns.  However  the  inversion  procedure  leads,  in  soae  cases,  to  dif¬ 
ficulties  because  of  the  presence  of  negative  exponents  tending  to  sero  in  the 
determinant  of  the  denominator.  Other  prograas  have  tried  to  avoid  the  inver¬ 
sion  problea  by  using  a  trial.  The  procedure  involves  selecting  values  for  B 
and  0.  and  solving  the  system  of  (An  -  !)  equations.  The  solution  is 

n 

evaluated  on  how  well  the  surface  conditions  are  aet.  A  second  pair  of  values 
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for  B  and  D  la  selected  and  a  new  solution  obtained  procedure.  Since  the 
n  n 

procesa  is  linear,  a  good  estinate  of  Bq  and  Dr  can  be  nade  by  linear  interpo¬ 
lation.  The  difficulty  lies  in  the  appropriate  choice  of  the  values  of  B  and 

n 

I>n  to  enaure  a  numerically  correct  Interpolation. 


However,  neither  of  the  approaches  are  appropriate  for  all  cases  of  a 
frictionless  interface.  The  vertical  def lectin  at  the  surface  of  a  two  layer, 
the  thickness  of  the  first  layer  being  H,  is  given  by 


I+u.  J  (mr).J  (ma)  . 

1  rm  o  1  , ,  2  -mH  „  2  aH 

p*“fT"  Jo  - 12 -  *A,B  c  "  *«■  « 


r 


’1“ 


-  [ 2-4u j+»H)  C^me”*H  -  ( 2-4p ^-mH)  da 


In  the  case  of  a  frictionless  Interface,  this  relation  becomes  (BURMISTER, 
19*5) 


-  pa- 


2( l~Uj )  .  ^(mrJ.JjCma) 


/; 


.  re2aH  -  (2F-l-2mH)  -  (l-Fle'™ 

F#2-1  ♦  (2F-1 )2mH  -  (1+2«2H2)  ♦  (1-P)e"2-H 


1  dm 


where 


(l-U^)  ♦  n(l-pj)  E2  (1  +  Wj) 

F  -  - —j- - - -  and  n  •  r— .  — , — - - 

2(1-w2)  1  +  v2 

For  ■  *  0  the  term  in  brackets  becomes  indefinite  (0/0).  This  condition  is 
not  slngif leant  when  computing  stresses,  because  the  product  of  the  Bessel 
functions  is  also  sero  at  the  origin  (JQ(mr).J j(ma)  *  0  for  m  »  0.)  However, 
in  the  case  of  deflection  the 


V 


a 

7 


llm 

m*0 
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As  a  result  it  is  necessary  to  determine  the  tern  in  brackets  for  n  ■  0  and 
thus  to  have  it  in  a  sufficiently  closed  font  to  be  able  to  compute  it. 

1.3.2.  Over  and  underflow  problena 

In  the  numerical  integration  procedure  m  varies  from  o  to  a  value  high 
enough  to  ensure  convergency.  Stated  differently,  the  integration  procedure 
can  be  stopped  when  the  terms  of  the  series  become  small  enough  so  as  not  to 
have  an  Influence  on  the  final  result.  The  value  of  m  to  achieve  convergence 
may  be  as  high  as  20  to  30.  The  nature  of  convergency  can  be  examined  by  con¬ 
sidering  the  two-layer  problem  developed  in  the  preceeding  paragraph. 

The  parameters  and  ,  from  which  the  values  of  all  the  other  parame¬ 
ters  can  be  deduced,  are  given  by: 

c _ [(l-F+mH)el,H  -  (l-F)e~"H) _ 

1  Fe2-H  +  (2F-l).2mH  -  (l+2m2H2)  +  ( 1-F)e~2-H 

D _ [Fe*H  -  (F-mme"**1] _ 

1  Fe2aH  +  ( 2F- 1 ) . 2mH  -  (l+2m2H2)  +  (1-F)e~2-H 

These  expressions  can  be  examined  in  terms  of  H/a,  where  a  is  the  radius  of 

the  circular  loaded  area.  At  an  H/a  value  of  5,  a  memory  overflow  will  occur 

aH/a 

for  values  of  m  above  10  as  a  result  of  such  exponential  terms  such  as  e 

and  e2-H^*.  However,  this  problem  can  be  overcome  by  dividing  both  numerator 

2aH 

and  denominator  by  e  . 

Expressing  the  parameters  in  this  form  will  result  in  an  underflow.  However, 
most  computers  have  a  routine  that  sets  variables  subjectsd  to  underflow  equal 
to  saro.  If  such  a  routine  does  not  exist,  it  can  be  build  into  the  program. 


-  16  - 


By  using  the  transforaed  relations  for  and  convergency  occurs  quickly 
and  in  a  coapletely  safe  way  (l.e.  the  nuaerators  both  tend  to  zero,  while  the 
denoalnator  tends  to  a  constant  F). 


This  can  be  obtained  autoaatlcally  for  the  two  layer  systea  by  writing  the 
boundary  condltlona  at  the  surface  (z  -  -H)  as  follows: 


A^e  +  B^e-"**  -  C ^ <  1— 2m ^  +  aH)e  +  D^( 


-2aH 

e 


1.3.3.  Stresses  and  Dlsplaceaents  at  the  Surface 


Convergency  is  slow  when  surface  stresses  and  dlsplaceaents  are  computed 
because  the  paraaeters  B ^  and  content  at  this  level  a  constant  term  in 
their  nuaerators. 


We  use  the  surfsce  conditions  to  express  B^  and  Dj  in  function  of  Aj  and  C ^ . 

Each  relation  is  then  split  into  two  parts:  an  integral  of  a  product  of 
Bessel  functions  which  is  coaputed  analytically  and  an  integral  containing 
only  negative  exponents  in  the  nuaerator  which  ensure  noraal  convergency.  In 
the  case  of  a  unifora  distributed  load  the  analytical  Integrals  are: 


J0(«>  (as )  da  -  ^  for  r  <  a 


■  u  f"  r  ■  • 

■  o  for  r  >  a 


for  r  <  a 
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“  2^  r  <  a 


Jj(ma) 

/  J  («r)  “  dm  which  is  a  special  case 


•  Ji(aa) 

For  r  »  o,  one  has  /  - - dm  -  1 


'  o  m 


„  JQ(mr).J  (ma) 

For  r  -  a,  one  has  /  — - - - dm  ■  2/w 


For  r  <  a,  one  has  J 


^  ^(mrJ.J^ma) 


dm  -  F  (1/2,  -  1/2; 


where  F  is  the  hypergeometrlc  function  of  GAUSS: 
-  (a)  (b) 

F  <*>  bi  Ci  *>  *  '  TT(c)  °  •* 
o  n 

(a)n  -  a  (a+1)  (a+2)  ...  (a+n-1) 


(a)o  -  1 


For 


J  (mr).J  (ma) 

■■  >*  i:- - r-5 - d. 


2r  F  (1/2,  1/2;  2;  £_) 


In  the  case  of  a  rigid  load  the  analytical  integrals  are: 


J"  JQ(«)  ain  (ma)  dm  -  (a2-r  V1/2  for  r  <  a 

•  o  for  r  >  a 

£  J  („,  .  t  ,.-(.W»l/2  for  r  < 


■  ~  for  r  > 
r 


fD  •*!<■*>  dm  -  £  [a-(a2-r2)l/2l  for  r  <  a 


■  ~  for  r  >  a 


V«>  ~ dm  -  w/2  for  r  <  a 


.  .  2  2-1/2, 

»  arctg  (a  (r  -a  )  J 
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1.3.4  Accuracy  In  the  First  Layer 

As  with  computation  of  stresses  at  the  surface  the  numerical  computation 
of  stresses  In  the  first  layer  or  near  the  surface  converges  slowly. 

Again  we  express  the  parameters  and  ,  associated  with  small  exponentials , 
in  function  of  the  parameters  A^  and  and  split  the  obtained  relations  In  an 
integral  of  a  product  of  Bessel  functions  and  an  exponential  (known  as  an 
LIPSCH1TZ-HANKEL  Integral)  which  has  to  be  treated  In  a  particular  way  and  an 
integral  with  normal  numerical  convergency. 

Infinite  LIPSCHITZ-HANKEL  integrals  can  be  transformed  In  finite  and  thus 
lntegrable,  integrals  when  the  Bessel  Indices  are  identical. 

In  the  case  of  a  uniform  distributed  load,  a  general  solution  can  be  achieved 
by  considering  the  stresses  and  displacements  in  a  semi-infinite  body  submit¬ 
ted  to  a  force  P. 

2 

With  P  «  pwa 


lim  pa  /“  Jo(mr)  J^ma)  dm  -  ~  JQ(mr)  dm 

S"0 


Thus  we  may  write  for  example 


p*  £  Vr)  **■  “  U  /o'  ^o  r.  P  V*>  *"**  dpde  d" 


2  2 

where  p  •  (a  +r  -2ar  cos  8 ) 


19  - 


In  this  case 


.  -mz  22-1/2 

J0  J0(np)  e  dm  «  (p  +z  ) 

so  that  the  Infinite  Integral 


pa  f”  J  (mr)  J.(ma)  e ~mz  dm 
‘  o  o  l 


Is  transformed  In  a  finite  integral 

p  ,2i  ,a  p  dp  d9 

2w  Jo  Jo  ,  2  2. 1/2 

(p  +z  ) 

One  of  the  two  integrals  can  be  solved  analytically.  The  required  integrals 
are: 


pa  J“  JQ(mr)  JjOm)  e-”2  dm  ■  2^  ft 


,  ,  2  V^2 

2z(a  -x  )  dx 


2*  J -a  /  22  2  .  .  ,  22  2  ,  .1/2 

(z  +x  +r  -2xr)  (z  +r  +a  -2xr) 


pa  /"  JQ(mr)  J^ma)  mze_,nz  dm  -  /*a 


2z(a2-x¥/2  dx 

(z2+x2+r2-2xr)  (z2+r2+a2-2xr) 1/2 


+  ^/a 
2w  J  -a 


2z  (a2-x2)1/2  dx 
(z2+r2+a2-2xr)3/2 


_  /  2  2  „  .  ,  2  2.1/2  . 

p  .a  2z(r  +x  -2rx)  (a  -x  )  dx 

"  2w  1 -a  ,  2.2  2  ,  ,  ,  2.2.  2  ,  T37I 

<z  +x  -*r  -2xr)  (z  +r  +a  -2xr) 


.  /  2^_2  ,  .  ,  2  2.1/2  . 
p  ft  4z(r  +x  -2rx)  (a  -x  )  dx 

"  2»  J-a  ,  2  2  2  ,  .2  /2.  2A  2  ,  TT72 

(z  +x  +r  -2xr)  (z  +r  +a  -2xr) 
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J  (mr)  J  (ma) 


_  - - -  -  .2  2  2  1/2  2  2  1/2 

_ .  r*  o  1  _-mz  j_  _  p  ra  (r  -2xr+a  +z  )  +  (a  -x  ) 

pa  - - -  e  dm  ^  /_a  In  -2-  -  -g  vl/2 - -  2  072  dx 


72^7  27  27172  :  2  2.1/2 

(r  -2xr+a  +z  )  -  (a  -x  ) 


and  following  LIEBSCHITZ-HANKEL  integrals 


J  (mr)  J. (ma) 
f“  1  1  -mz  . 

pa  J  -  e  dm 

J  o  mr 


a  fii  1  ,  z 

**7“  Jo  2  1  2  2  1/2 

w  (z  +w  ) 


sinOdO 


2  2  2 

with  w  **  a  +r  -2arcos  0 


J  (mr)  J  (ma)  2 

r<*  1  1  -mz  .  a  fir  1 

pa  jQ  - - -  e  dm  -  p—  J  T  "  2"! 

(z  +w  ) 


sin  9 d6 


2  2. 

pa  /“  Jt(mr)  J^ma)  raze-"12  dm  -  3p^~-  — j — 2~5/2 

(z  +w  ) 


sin20d0 


In  the  case  of  a  rigid  load  we  expand  either  the  Bessel  function,  J^imr) 
or  Jj(mr),  either  the  sine  function,  sin  (ma),  depending  on  the  relative 
values  of  r  and  a,  and  solve  the  resulting  series  of  integrals. 


The  required  Integrals  are: 


J  (mr)  sin(ma)  e 

/**  -2 -  dm  - 

'o  m 


if  r  >  a 


E  TO TT  T’ZITTZnVl  F  n+1/2’_n5  U 

o  fa  *z  )  r  +z 
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if 


r  <  a 


arc  tg  (a/z)  +  az 


“  (-l)n  r2n  r  (n+1/2) 
1  T(l/2)n!  (a2+z2)1/2 


2 

n+1  ,-n+l ;  3/2; 

a  +z 


j  J  (mr)  sin(ma)  e  mZ  dm  = 
'  o  o 


if  r  > 


•  ,  isn  2n+l 

°  ,  2  2.  2 

(z  +r  ) 


2n+3  , 

~2 — »  “*»»  l; 


2  2 
r  +z 


if  r  <  a  2a  £ 


(-l)n  r2"  r  (n+3/2) 
o  r(l/2)  nf  (a2+z2)n+1 


n+1,  -n;  3/2; 


2 ,  2 
a  +z 


/  J  (mr)  sin(ma)  mz  e  dm 
'  o  o 


if  r  >  a  az 


0  2  2 
2z  -r 


(-l)n  a2n+1(n+l) 


_  _  *  4.  0_  y  \  *  /  O _ 

(:V)5«  1  .  ,  3S1 


(z2+r2)  2 


2n+3  , .  , .  r 

2  »  ~n-i,  i,  2  2 

r  +z 


Ur«  r2n<"t‘>  >  , 

o  n!  r( 1/2)  (z2+a  )n  2 


n+2,  -n;  3/2; 


2  2 
a  +z 


-mz 


J. (mr)  sin(ma)  e 

r~ - dm 

•  o  m 


if 


r  >  a  -  - 
r 


•  ,  ,xn  2n+l 

az  .  rz  (-1)  a  _ 

rV+z2)1'2  2  \  _ 

(r  +zZ) 


12+3  _n+1.  2. 

2  »  l,Tl  >  2 


r  +z 


if 


r  <  a  a 


•  /  ,  \n  2n+l  _  . 

£  — Izll  £  r  (n+3/2) _  p 

o  (n+1)!  T(l/2)  (z2+a2)n+1 


n+1 , 


3/2; 


2 

a 


2 

a  +z 


2 


J^Gnr)  sin(ma)  e-02  dm  - 
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»  ,  ,.n  ,  2n+l  .  2 

r  v  „  r  (-D  (n+1)  a  .  r  2n-(-3  r 

if  r  >  a  S  2^3  F  2  ’  n;  2’  1~2 

o  2  2  2  ** 

(r  +z  ) 


if  r  <  a  21  ^~1)-  - - aZ  2*1+2 }  F  n+2*  -n;  3/2;  "TT 

o  n!  T( 1/2)  (a  +z  )  L  a  +z 


f  J.(mr)  sin(ma)  mz  e  mZ  dm  = 
1  o  i 


r  v  n  .  r  (-l)n  a2n+1.r.z  (2n+3)  (n+1)  2n+5  .  r2 

if  r  >  a  z  £  - 2n+5 -  F  ~T~'  'n;  2;  ~2~2 

0  2  2  1“  r  +z 
(r  +z  ) 


if  r  <  a  —  23Z2~3  }  +  2raz  E  - - —  F  n+2,  -n-1;  3/2; 

(a2+zV  1  n!  T( 1/2)  (a2+zV  a2+z2 


PART  2:  NUMERICAL  RESOLUTION  OF  A  FOUR  LAYER  SYSTEM 


2.1  THE  MATHEMATICAL  ANALYSIS 

As  shown  in  art  1.3.1,  the  mathematical  analysis  must  be  different 
depending  on  the  layer  interface  conditions. 

2 . 1 . 1  Full  or  Partial  Friction  at  the  Interfaces 

The  case  to  be  solved  for  first  is  full  or  partial  friction  at  the  inter¬ 
faces.  Boundary  conditions  for  this  case  will  be 


-  At  the  surface 


rzi  +  1 
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ui  =  Xiui+1 

where  is  a  factor  for  partial  friction. 

-At  the  bottom  w  =  0 

or  =  0  if  rigid  bottom 
A^  =  =  0  if  4  -  th  layer  is  semi -infinite 

Consider  an  isotropic  four  layer  problem  which  will  result  in  a  system  of 
16  equations  with  16  unknowns  (A^,  D^).  This  system  of  equations 

is  difficult  to  solve  analytically  in  the  same  way  that  BURMISTER  did  for  a 
two  and  three  layered  system.  Using  the  same  approach  results  in  mathematical 
errors  introduced  in  eliminating  unknowns.  As  a  result,  it  is  necessary  to 
modify  the  analysis  procedure  to  solve  accurately  the  numerical  problems 
detailed  in  Part  1.  The  main  objective  of  the  approach  is  to  obtain  an 
expression  for  each  unknown  parameter  consisting  of  a  numerator  containing 
negative  exponents  only  and  a  denominator  containing  a  constant  term  and  nega¬ 
tive  exponents.  During  the  integration  procedure,  when  the  variable  tends  to 
infinity,  the  numerator  will  then  tend  to  zero  and  the  denominator  to  a  con¬ 
stant  value.  In  the  analysis  the  exponents  must  appear  in  close  form.  The 
factors  multiplying  the  exponents  may  then  be  expressed  in  a  more  comprehen¬ 
sive  form.  Sequential  steps  in  the  mathematical  analysis  for  a  rigid  bottom 
are  as  follows. 

First  Step 

Replace  in  the  boundary  equations  of  the  third  interface  the  parameters 
A^  and  by  their  values  obtained  from  the  rigid  bottom  condition. 
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Second  Step 

Write  all  the  interface  conditions  in  matrix  form. 

T  T 

-  At  the  surface  MI  (A^  D^)  =  (I  0) 

where  Ml  is  a  2  x  4  matrix 

T  T 

-  At  the  first  interface  M^(A^  Bi  ^1  =  M2^A2  B2  ^2  U2^ 

where  M^  and  are  4x4  matrices 

T  ,  T 

-  At  the  second  interface  M^(A^  D2>  =  M4^A3  B3  ^3  ^3) 

where  M_  and  M.  are  4x4  matrices 
3  4 

T  T 

-  At  the  third  interface  M5^A3  B3  C3  °3^  =  M6^B4  U4^ 

where  M,  is  a  4  *  4 
and  is  a  4  x  2  matrix. 

Third  Step 

Invert  the  matrices  and  M^. 

The  system  becomes 

MI  (Ax  Bt  Cl  =  (1  0)T 

U,  »,  Qj  d/  -  m;1  .M2(i2  B2  C2  D2)T 

(a2  b2  c2  d2)t  =  rt"1  .m4(a3  b3  c3  d3)t 
<*3  B3  C3  U3)T  *  “I'-V6*  D4,T 
and  finally 

mi.m~1m2.m”1.ma.m~1.m6(b4  d4)T  *  (1  o)T 

The  product  MI.M^,M2.M3*.M4.M3*  .M^  is  a  2  x  2  matrix,  so  that  we  can  write 
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and  solve  the  resulting  system: 


=  _ *12 _ 

4  all *a22  “  a21 *a12 


4  a11,a22  "  a21 *a12 


Utilizing  the  different  matrix  equations,  we  express  then  the  other 
parameters  as  functions  of  and  D^.  To  obtain  the  exponents  in  close  form, 
the  matrices  are  split  in  such  a  way  that  the  exponents  can  be  factored  out  of 
the  brackets.  For  example,  in  the  isotropic  case: 


M 


-1 

1 

M2 

-1 


M 


M4 

-1 

5 

M, 


"  4(1-Wl)  [  e  ,M11  +  6  *M12 

[eX.M2l  +  e  X.M22] 

_  4(1-u2)  [  e  ,M31  +  6  *M32 
[  ey.M41  +  e'y.M42  ] 

“  4(l-p3)  [  e  *M51  +  ®  ,M52 

[  e  .M61  +  e  .M62  ] 


where 


mh. 


y  *  mihj  +  h2) 
z  *  m(h1  +  h2  +  h3) 
u  “  m(h^  +  h2  +  h3  +  2h4) 
h^  being  the  thickness  of  layer  i. 


] 


] 


J 


The  terms  of  all  the  matrices  are  in  closed  form  and  can  be  introduced  as 
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k 


\ 

r 


inpuC  values  for  the  numerical  procedure.  The  terms  of  the  resulting  products 
are  not  expressed  in  closed  form,  however  it  is  sufficient  to  know  which 
columns  or  rows  of  each  matrix  contain  only  zeros.  In  the  final  product 
intermediate  matrices  disappear  because  they  are  identical  zero.  As  a  result, 
those  matrices  preceeded  by  positive  exponents  vanish.  Subsequently,  the 
values  of  all  the  unknown  parameters  except  for  the  parameters  Bj  and  Uj  at 
the  surface,  can  be  expressed  in  terms  of  a  numerator  with  only  negative 
exponents  and  a  denominator  with  a  constant  term  and  negative  exponents. 

At  the  surface,  the  values  of  the  parameters  B^  and  contain  a  constant 
term  in  both  the  numerator  and  the  denominator  followed  by  negative  exponents 
which  precludes  obtaining  satisfactory  convergency.  As  an  alternative,  the 
parameters  B^  and  0^  are  expressed,  utilizing  the  surface  conditions,  as  a 
function  of  the  parameters  A^  and  .  The  resulting  expression  for  a  given 
stress  or  displacement  becomes: 

°  ■  pa  JQ  J^mO.J^ma)  l  tt  +  f^A^)  +  f^Cj)  J  am 

where  K  is  a  constant  term  and  f^(A^)  and  f^(C^)  are  functions  of  the  parame¬ 
ters  and  which  converge  normally.  Splitting  the  above  relation  into  two 

parts  results  in: 

o.  *  pa  j  J„(mr) .J, (ma) .  K  dm 

1  '  O  U  1 

which  can  be  solved  analytically 

°2  “  pa  /q  Vmr)*Jl(ma)  l  fl^Al^  +  f2^Cl^  ^  dm 

which  can  be  solved  numerically. 


L 
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Near  the  surface  the  previous  probJ.ee  of  and  converging  slowly 
arises  too.  Therefore,  utilizing  the  surface  conditions,  the  expressions  for 
stresses  and  displacements  are  split  Into  two  parts.  The  first  part,  known  as 
a  LIPSCHITZ-HANKEL  integral  (paragraph  1.3.4),  is  solved  analytically  and  the 
second  part  la  solved  numerically. 

2.1.2  Full  Slip  at  the  First  Two  Interfaces  and  Full  Friction  at  the  Third 
Interfaces 

A  different  procedure  must  be  used  to  mathematically  solve  a  layered  sys¬ 
tem  with  frictionless  conditions  at  some  Interfaces.  The  procedure  described 
in  previous  paragraphs  cannot  be  applied  here.  The  boundary  conditions  for 
the  case  of  full  slip  for  the  four  layer  system  are: 


U) 

(2) 

(3) 

(4) 


o  for  a  rigid  bottom 
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The  systeaa  of  equation?  at  the  first  two  interfaces  cannot  be  expressed  in 
aatrix  fora  because  two  of  the  four  equations  are  hoaogeneous. 

Sequential  steps  in  the  aatheaatical  analysis  are  as  follows: 


Firat  Step 

Replace  in  the  boundary  equations  of  the  third  interface,  the  paraaeters 

A  and  C  by  their  values  obtained  froa  the  rigid  bottoa  condition. 

4  4 

Second  Step 

Write  the  boundary  equations  at  the  third  interface  (full  friction)  in 
aatrix  fora 

T  -1  ,  T 

(Aj  Cj  Dp  *  *M^(B^  D^) 

Third  Step 

Using  the  surface  conditions,  express  A^  and  B^  as  a  function  of  and 


(*,  b,)T  -  Ho(C,  U,)T 

Using  conditions  (l)  and  (2)  at  the  first  interface  replace  and  Bj  by  their 
values  and  solve  the  systea  by  expressing  and  U^  as  a  function  of  A^, 

^2*  ^2*  ®2* 


(Cj  d,)t  .  H,(*2  b2  C2  D2)T 


Using  condition  (4)  at  the  first  interface,  replace  A^,  B^,  and  by  their 
values  expressed  as  functions  of  A^,  B^,  C^.  and  Using  conditions  (3)  and 
(4)  at  the  first  interface,  express  A^  and  B^  as  a  function  of  and 
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(A2  Bj)  ■  ^2^2  °2^ 

Using  condition  (l)  and  (2)  at  the  second  interface,  replace  and  by 
their  values  and  solve  the  systea  by  expressing  C  ,  and  aa  functions  of 

A3,  C3,  ^3* 

(C2  D2,T  *  M3<*3  '3  C3  D3)T 

Using  condition  (4)  at  the  second  Interface,  replace  A^,  B^ ,  ,  and  by 

their  values  expressed  as  functions  of  A^,  B^,  and  D^.  Conditions  (3)  and 
(4)  at  the  second  Interface  are  reduced  to  the  following  systea 

M4(A3  B3  C3  D3)T  -  (K  o)T 

where  is  a  2  x  4  matrix  and  K  a  function  of  the  integration  variable. 
Finally, 

M4.m"1.M6(Ba  VT  -  (K  o)T 

a  system  which  can  be  solved  as  in  previous  paragraph. 

The  same  procedure  is  also  utilized  to  express  all  the  unknown  parameters  as 

functions  of  B,  and  D,  . 

4  4 

At  this  point  a  supplementary  difficulty  for  the  full  slip  case,  regard¬ 
ing  the  expression  of  the  vertical  deflection  must  be  considered.  The  verti¬ 
cal  deflection  is  given  by: 


1+u.  J  (mr).Jt(ma) 

w  -  pa  X1  fo  - - ST -  fi(AiBiciDi>  d" 


For  m  •  o 
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lia 

a*o 


J  (ar).J  (aa) 
o  1 


2 


and 


11a 

a«*o 


fi(AlBlClDl) 


0 

0 


which  result  in  the  above  expression  for  the  vertical  deflection  being  unde¬ 
fined.  However,  at  the  rigid  bottoa  w«0  for  all  values  of  the  integrating 
paraaeter  a^ ,  including  s*o.  Consequently, 

At  the  bottoa  f.(A,  B  C  D  )  *  o 
4  4  4  4  4 

At  the  third  Interface  f  (A  B  C  D  )  -  f  (A  B  C  D  ) 

33333  44444 

At  the  second  interface  f2^A2  B2  C2  ”  f3^A3  *3  C3  °3^ 
and,  at  the  first  Interface  f^(Aj  Ci  D1 )  “  ^2^2  B2  C2  D2^ 

Thus,  in  general  for  all  layers  f^(A^  ®i  Ci  V  "  °  ^or  *  **  °* 


Of  course  this  is  not  any  aore  true  in  the  case  of  the  last  layer  being  a 
seal-infinite  body.  At  the  origin  of  integration  (a"0)  we  then  write  the 
boundary  conditions  as  follows: 


-  1 


B5  +  Dj  (1  -  2u5)  -  1 
0  -B5  ♦  2u5  D5  -  0 


We  obtain  B^  ■  2 


- 


and 


1 
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so  that  we  find  the  origin  tern  of  each  deflection  by 

1+M. 

l(At  -  )  -  2(1  -  2Ml)  (^  -  D^]  -  l-b5  ♦  2  0-2u5>  »>5 1 


.2(1  -  n5) 


2.2  THE  NUMERICAL  PROCEDURE 


In  general,  nuaerlcal  Integration  can  be  accomplished  using  Simpson's 
method  of  the  form: 


l  U.b)  -  3  U„  ♦  «,  *  U2  ♦  ...  •  «2„.,  ♦  <2„l 


where  la,b)  is  the  Interval  of  integration  of  1  (a,b)  subdivided  in  2n  equal 
segments  of  length  h. 


However,  stresses  and  diaplaceawnts  are  obtained  by  solution  ot  expressions 
such  as : 


o  -  pa  /  J  (ar).J.(aa)  f(m)  da 
'oo  l 

where  the  interval  of  Integration  goes  until  infinity.  The  nuaerlcal  computa¬ 
tion  is  Interrupted  when  the  value  of  the  function  f  (a,z)  (in  fact  the  values 
of  all  the  parameter  becomes  smaller  than  the  imposed  convergency 

level  (i.e.  ^(m.r)  <  e).  The  rate  of  change  of  the  function,  ^(sk),  Is  sig¬ 
nificant  for  small  values  of  the  integrating  parameters  a  and  less  significant 
for  high  values.  A  reasonable  approach  appears  to  be  to  increase  the  incre¬ 
ment  of  a  (i.e.,  h)  for  higher  values  for  m.  Using  this  technique,  the  above 


relation  can  then  be  written: 


-  32  - 


I  Io,»]  -  yi-  IfQ  +  4f1  +  2f  J  +  f^]  for  m  <  Lj 
h2 

♦  y~  +  4f  ^  +  2f^  +  4fy  +  fg]  for  <  ■  <  Lj 
h3 

+  3~  ^8  +  ^f9  +  2f  10  +  +  *12^  for  L2  ^  “ 

For  practical  reasons  we  take 


1  |Li-  1W  miT  “fi..  *  2fi«  *  -  ‘Wi  *  2fi«„' 


The  length  of  the  integration  segment  can  be  modified  further  with  a  value  of 

A  f 

T"  ’‘i+in 


The  initial  length  h  of  the  integration  segment  is  chosen  by  the  user  and 
is  multiplied  by  a  factor  of  2  when  the  values  of  f(mz)  become  smaller  than 
10  ,10  ,  10  and  10  .  The  final  convergency  level  is  10  **. 

9 

For  the  solution  of  vertical  deflection,  fQ  *  0  and  is  computed  in  an 
appropriate  subroutine. 


2.3  COMPUTER  PROGRAM 
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A  computer  program  has  been  prepared  that  runs  efficiently  on  a  personal 
computer.  The  main  steps  of  the  prograa  can  he  resumed  as  follows. 

2.3.1  Prograa  Input  Procedures 

“  Input  of  the  data  (by  display  or  file)  for  loads  and  pavement  struc- 
ture 

-  Input  of  the  coordinates  (by  display  or  file)  of  the  points  where 
stresses  and  displacements  are  to  be  computed. 

-  Choice  of  the  length  of  the  integration  segment. 

2.3.2  Computation  Features 

-  Vertical  deflection  at  the  surface  and  in  the  first  layer,  for  m  - 
o. 


Parameters  B^  C^  for  each  value  of  m. 

Bessel  functions  for  each  value  of  m. 

Stresses  and  displacements  in  cylindrical  coordinates  for  each  value 
of  a. 

Stresses  and  displacements  in  cartesian  coordinates  for  each  value 
of  m. 

Vectors  containing  the  results  using  Simpson's  rule. 

Convergency  test  for  each  f  function. 
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-  When  the  final  convergency  of  the  numerical  part  is  reached,  compu¬ 
tation  of  the  analytical  part  of  the  solution  for  the  surface  and 
the  first  layer  must  be  added  (LIPSCHITZ-HANKEL  Integrals). 


2.3.3  Output  Procedures 


stresses  a  a  a  t  t 

x  y  z  xy  xz 


T 

yz 


displacements  u  v  w 
x  y  z 


-  principal  stresses 

principal  strains  z^ 

linear  strains  e  e  e 
x  y  z 

2.3.4  Capabilities  of  the  Program 

The  number  of  circular  loads,  with  different  radii  and  contact  pressures, 
is  limited  to  20  and  stresses  and  displacements  can  be  computed  at  30  loca¬ 
tions  in  the  horizontal  plane  and  at  each  place  at  30  depths  (included  8 
values  at  the  surface  and  interfaces). 

2.3.5  Software  and  Hardware  Requirements 

The  program  is  written  in  FORTRAN  77  and  runs  on  IBM  PC's  equipped  with  a 
8087  Math  Comprocessor  and  compatible  equipment.  The  executable  version 
requires  200  kilobytes  of  memory. 
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ALGEBRAICAL  ANALYSIS  OF  I90T10PIC  LAYERED  SYSTEMS  WITH  FIXED  BOTTOM 
AMD  PARTIAL  FRICTIOM  OOHDITIOHS  AT  THE  INTERFACES 


APPENDIX  1 


ALGEBRAICAL  ANALYSIS  OF  ISOTROPIC  LAYERED  SYSTEMS  WITH  FIXED  BOTTOM 
AND  PARTIAL  FRICTION  CONDITIONS  AT  THE  INTERFACES 

SYMBOLS 
We  write 

A-,.  A .  *\v»x  C;  -- 

Fw *  Fw=  A, 

ww ,  v<w  -  [^H 

E3  01+>k*)  Et> 

lw  *  *  A* 

E  A  (*+  hO  Ea[»*K-0 

X  =  >v>  H , 

y  «  Yv»  (^i  +  Ha') 

Zcw>(U,vHa  +  H3J 

w  (Hft-  Hi  +  +H^3 

where  Hj,  Hg,  H^  are  the  thicknesses  of  the  successive  layers. 
The  index  1  applies  to  the  surface  layer. 

We  shall  successively  analyse  a  two  layered,  a  three  layered  and 
a  four  layered  structure. 


A  1.2 


Chapter  1.  The  two  layered  structure. 

1.  The  boundary  conditions. 

Boundary  conditions  at  the  surface  (2  =0): 

<Tt=  A,  ♦  £>,  -  C,  (>1-2^0  +  35,  -l 

A,  -  £>,  ■*  C,  •  ■r  35,  2.^/.,  s.  \a 

Boundary  conditions  at  the  interface  (2  =  H^); 

Tz  A,fcx  r  2>,  e  x  -  C,  (>-  2f*.  -*  )  e.*  +  35,  O  V-  ♦  0  «•"  *= 
A2<->  *  b2  e-*  -  c4  (/t-2>^*-,0*X  4  £>2(4-  + 

tV*-.  A,«*  -2>,  **  ^  cr» +■  ^  C  2-^,-^)  *  X  *■ 
Ai**  -  r  Cl  *) «.*  +  D*  OM"*)1 
w:  a,**- 

u  •.  A,*’*  b.C*  +  C,(-^x)^. 

F^^Aie*+  ***'  *c*omK- 
Boundary  conditions  at  the  bottom  (2  =  Hj  +  H2): 

W  •.  AacV_  ^y_  CiO-Aft-y)^- 
c*«  o 

2.  Resolution  of  the  system  of  6  boundary  equations. 

We  write  the  conditions  at  the  interface  in  matrixform 

and  invert  matrix  Ml 

(A,  a,  c,  t>,JT=  Mi'.Mz 
after  having  replaced  A2  by  its  value  at  the  bottom 


A^=  Sh 


*  «- 


■rUY-*) 


A 


A  1.3 


Vf  A  -X 

n,  a - - 

-4^-hO 


A(a-  }*..) 


TV* 


i  (Mh+T) 

-(*/“*)  i  (l-^ivf) 

t 

Fw 

Tu.  ?v.(i-A^y) 


~t 

*  e 


(A-lyx.+*) 

{xy*->-) 

-Tw{r-Ayi*y) 


(A,6,C,S,)T»  -  J - [»*M„  ♦  «.***][  DOT 


a  .  J__  \ *%„!,  * w  „22  +  tzb  x;n)i2l  *-  Maiil  (Bt  pot 


4^-vO 


4-  * 


,  o  o 


Mu2JL  a 


4  ♦ 


0  o 


1-  V 


0  o 


f 


0  o 


+  + 


o  o 


+-  *- 


A  1.4 


«e  —  t>  +•  >4  y* ,  T\*.  ~  2^«.4  ^  T  w  -  Tv.} 

^4  *  *  1  Fvy  >  Fw .  |m  -  ^Ia. 

?i  =  Fw-K 

«  -  >\  +-  t  -  S)»-»yi  -  F»*»  »  S^'y^tFvv  ♦  Tv  -  2^vTv. 


_A _ 

A 


c  Mnat  O' ')  +■  ^  »2x  (•> 0^  &i*> 

*  £  ^ w at  (*t2.)  +  F)  naa  (j**)  J  Pi%.^  |j 


__A _ 

M^-hO 


—  i  J  ^ 

A(*hO  [  J 

4  \  Fiijj  Di  ^ 

-  "  l 

Ma*m)  j 

+■  ft*.  ♦  [_Rxi.  ■♦■  Di'j 

c,a-  4  ) 

A  h-h')  ( 

[*wx,  O')  e  ^  ■♦■  )l«»x  C**')J  bfZ^* 

+  [M,i,C>.*K2(,'~y)  *  Hnll  (>,t)]  j>2i* 

-  4  [ 
A  (i-V..)  L 

^1*  "**  4- 

1\  _  .  A 

*M  -*  -  »■  , 

4(4-H 

^  M  ioai  O')  *■  4  FaiJ 

4  (A *•")*■  ^  f  3>» 

>1 

/ 

j  bi  4  ^  fUz.  4  D?  ^ 

(A,  lb,  C,  J,)T»  -  -±  —  MS  t  MFJ  ®z)T 

We  write  the  boundary  conditions  at  the  surface  in  matrix  form 
Mo  (A,  V9,  c,  J>0Te  (*  «)T 


-  \ 


s 


I 


-  (a-  2^. ,)  2$a>) 

2)*t  2y<, 


4 


f 


A  1.6 


m..  t*t.  t>oT*  (-*(->-hO  »)T 


A«  A,x  1 

Ff„ 

FF,x 

W..Hfta 

VI*  NF  = 

An  An 

FFx, 

FV»x 

A„» 

ft,,*  •*■  ftx^ 

*  l^O  1)1 1  *  ^-2^«i)  ft/»i 

A,x  • 

ftrt£^  4  ftzt 

—  (-4— 2f*t)  4  ftAz 

Ax> « 

ftn^k  -  flz, 

4  2^.,  £  k  4  2^.  ft^A, 

AjX  »•  ft,!*^  -  ftl2  *  2^*.,  ftw*1*  +  2J^i  ^Az 


?F„  *  rM  ♦  (4-2^0  *» 
fF,x  *  F»X  4*  («*-2j**)  Fax 
FPx,  »  -?z,  *  Fa, 

fftxi  r  -  Tzx  4-  2^-*  FAx 
We  develop  the  matrix  equation 

(A„  4-FFw)  &2  *  (  A,x  ^F,x)  3j  «  -  -4(*-F') 
(An  ♦>  F^xi)  bx  +  ( Axi  •♦■  Ffix)  Dx  «>  s> 
and  solve  the  system 


^5 2.  =  -  4  (m- h') 

T*r  *  A  C^-h') 


Axx  ??zx 

V 

Ax,  *  F^a) 
V 


V*  Am  (  Aja.  4  FVxz)  ~  A»x  (  Ax,  4  Ffz,)  4pf»  Ajt 
-?F>xAj,  *  ??„JF?xr  -  FP>x-Ffx, 


A  1.7 


We  develop  FF11.FF22  -  FF12.FF21 


FF*.FFW  -  FFj.1  -  F^t  * 

A1  the  other  members  of  the  denominator  contain  negative  exponents. 

The  limit  value  of  the  denominator,  for  m  =  *  ,  is  equal  to  a  constant 
F1.F2  +  F3.F4 

3.  Relations  for  the  parameters. 

The  values  of  the  parameters  and  Cj  are  deduced  from  the  matrix 
equation 

/A,  b,  c,  a,;Ta  -  \  m  +  MFA  Cbi  Di)T 

1  J 

(A,  o  C,  o  )  T  =.  — - -  HK  Cb,DOr 

a  (*-H 

A„^Azi  +  mx)**  _  Kn  (A»»  +  FFn)e> 

A  ,•  - 

V 

x  An  £Aau  +  ?F27.)  -  A-ix.  (An  *  ffzi)  «.  * 


(A**  +. TfzxJC5’'  - 

Ct  «  * - - 

v 

c  m  V (A 22  +  ??».)«:*-  »pe*Q«t 


The  values  of  the  parameters  Bj  and  are  deduced  from  the  surface 
conditions 

bi=  2^.,  ♦  A, 

1)|  ~  4-  2.A)  t-  C'i 

n,  -*  .  r  ,  .  v  Aaz.  +•  FFZ2  -X 

©it,  -  -4h-Ui)  -  e. 

7 

-/.[.-vl  »-* 

7 


.aft  _ . . I ... 


A  1.8 


Chapter  2.  The  three  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  second  interface  (z  =  +  H2): 

At*?  *  T>i 

W+  b>«rV-C*  (>,.2^j-y)e.V  v 
X«  •  ^t*V-  *>z«rY  ♦  C*C2-HM)eV  v  3)a  C2-K-v)e-7. 

Vv  :  Ate*-  b2e'v  -Cx(Z-4)n."y)  -  T^z  e 

k»v  I/***-  b»e'y  .  C*  £z  -Ay  -  D»0 ->$)*» 

:  Ai*y  V-  ftj  eTy  ^  ^i+y^ey  _  Da.  £-i-y )e~y  a 
kw  ^Aj*?  ♦-  ba«~y  +  + 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  N2  +  H3): 

W  :  *  -  C$  j  -%)  f  *0 

C  ^  =  o 

2.  Resolution  of  the  system  of  10  boundary  equations. 

We  write  the  conditions  at  the  second  interface  in  matrix  form 


(Ax  JbiCx  t>x)T.  J>j,y 

and  invert  matrix  M3 


(At  KCt  1>i)T*  D»)t 

after  having  replaced  A3  by  its  value  at  the  bottom 


A^**e 
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We  write  the  boundary  conditions  at  the  first  interface  in  matrix  form 

n,(A,  to,  C,  1>0T«  Oi*i  Ci- 
and  invert  matrix  Ml 
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Ml-1  is  given  in  §  1.2 
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We  write  the  boundary  conditions  at  the  surface  in  matrix  form 
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3.  Relations  for  the  parameters. 


These  relations  are  deduced  from  the  different  matrix  equations 
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Chapter  3.  The  four  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  third  interface  (z  =  ^  +  H^) 
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2.  Resolution  of  the  system  of  14  boundary  equations. 

We  write  the  conditions  at  the  third  interface  in  matrix  form 
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»e  write  the  boundary  conditions  at  the  second  interface  in  matrix  form 

M*  Ma  (A*  H)T 

and  ivert  matrix  M3 
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3.  Relations  for  the  parameters. 

These  relations  are  deduced  from  the  different  matrix  equations 

A»«\  -  -4—  .  p„-  ^e*  *  ?,x.  X>a£*1 

^  ^(.Px,  4-Lxi^)  bAt^  V  (  Fxi  ♦-  Lxx}  Da&^J 

- - - -  i  *■ 

_  _i — . "  {K  *  Ui)  Wv  *  Cp-S2  1-U*}  Rsf^ 

*G«-h)  ' 

Al^s  -  [w»  (*>')  '*  V'J  -V„(c».*)t"t',t'>’;r  t|*(J»t'0] 

CiJ.  -_i -  lv*(A,{)C1’’^ v 

X*'1* — —  I  vai (0*0 
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The  computation  of  the  parameters  B2e"x  and  D2e"x  needs  that  of  the 
intermediate  parameters  A3e2y,  B3X,  C3e2y  and  D3e~x: 

O.  -  - [  *.  ft,  v  D, 

V**  -  -4 -  V*  + 

-  — -- — .  J  (  p*,  4-  u.)  ft* *-*  *  Ljt,x)  i>a  €"*] 

*(h*0  L 

Bl** °  '  iiZTT)  I'®*'  (*>''')*'  •>•  *f«  te»")  <■ **» <?*&)*-"  *  vM  V"’J 

Dlt  «  -  — .  [dUi(A*e*,)e*  4-  kvjj.  ♦-  0a3  (Ci*' )«.”*  +■  JcVv  ^**1 

u  -1 

A*’*  L  F"  *a(®«0  -  s„  (w*>J 
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4.  Determination  of  the  stresses  and  the  displacements. 

4.1.  Mathematical  procedure. (in  the  case  of  a  flexible  load) 

The  program  calculates  the  stresses  and  the  displacements  using  following 
relations: 


_  Vw«. 


IT**  ^  I  ^Ai**1*1''**- 

o  ypfz.  ^ 

r«  _  _  Vi  I  y.  (w)  • 3'  (■“•“)  +-  ’b  .k-V 1 1,1 

"  »>  *  1  *” 


(Tr-f© 


4 1  *y*r  *•  y^X  1 


Xrzs-Va 


t  -Vhi 


+C;w.(l^  H'1)''*  S’>v’  ^  2)"  U'2  ‘l  J  * 


vv.^iH.yo. 


=  _>hH. 


E, 


■r 


*  T4»f)-5H  \  A\v*V"a-  W< 

j  **  *  4  N-k'Mj 

-Op-Ay.-**)*-  r  J 

^,C»r).?,H  j  AiwNTS 

'  ,  ywx  ^yv,^-^e  UK, 

<-Cstv»QA^*>  ~  V  J 


Stresses  and  displacements  are  calculated  in  a  system  of  cylindrical 
coordinates.  The  stresses  due  to  several  loads  are  to  be  added  together 
Therefore  we  must  express  them  in  cartesian  coordinates  using  following 
relations: 


Tx  -  <rr.  cci'oi  -f 

O' 0  .  tV%n"Ci  a 

Tr-*-  [%  ^ 

1 

7r-n 

* 

tv'  2.  k 

Oy  =  0~T-  y^«  +■ 

^  .  CfW  a. 

JL 

rr-  r* 

JL 

to  2* 

Tyx  -  XTR. 

a  T-r*  •  * 

"toof  *  (3>  -  ^ fe) 

t^TiOt-  /Vv»  £*  - 

5>-0V 

1 

/vU  Jl* 

lA*  s  b*-f.  M** 

Yy  «.  lm\  Oi 

W  ®  W 

wherein  the  signification  of  the  angle  «  is  illustrated  below. 
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4i2<1  _Stresses_and_dlsglacements_at_the_surface_£z_=_0) 


-  =  j> 

=  */x  Cr=et) 
-  O  Cr?^) 


1*0 

T,(v"*)lA,vb,+  (^^hOc'  ~ 

u 

Q 

»-)*•£-  3"o(*ir)  Y  (.4+2h0  4 ^ 

»*o 

_  (*►»*)  (>»+2^o<.i~ 

-?)>»  (-!*)-•)  X. l».r) I, +  2l.,(c,«.y)J^X*fc» 

with  („  1^)1  rr'0 

*  (  o  tr>) 

_  ,  .  >«:  3,K)[A,  ♦  t>.  ♦  c.  - 

Jo 

>.yw  V*  Ttl^O-  !)(**)  fc.*C\ -3),J  *U~ 

Jq 

_  la 

v  r,  2^*1* 

-t>^-y..)j  V“r>  - 


.4  ( 

La~1^')  I  ^  far)  •  3*  fa*)cLi*+*  \  ^1-2^,*)  £-  (_r-  c\) 

l  *  (0*0 


?,  (y»r) • 
Tnr 


-(4-2^.)l  C r^) 

(.f^ 


—  trx« 


o 
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E 


■'o 

V^V  y~  j  y(.»O  TiH  ^(A,<f)  *■ 


!*  T.^mO.T,  (»>«) 


i  x  Lr<") 

v 

>«  a.  (.<7<M 


4i3i_S  tresses^anddi  sgl  acerrents  _i  n_the_f  Irs  t_  1  a^er_£_0*^K_<^_Hi_2 ; 

-  fx.  jao.  T% (i^r) . T, (w> +)  -  c« 

-►i,  ^ 


jjtv  3,o»0.  *(*••)•  c^u*'  * ^ol^O ^,CY^,v)v''' e 


-(x-vAO 


+  >  JV«rV T.C-1 


4 

T»  (»v»r)  T,  (*»<*)  ^  * 


y  if  t  t»*~-  ^ 
*  V  **4* 


=  .t 


®  *T  rv,VAr'-2xrX>»t+f,+  e''  -Zxr  )* 


(/  =  «) 


(rtc>> 
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t* 


-WvV 


T«{r»r).  x  v,.  _»*_ 


(r~) 


.± 


i  y/» 


**  J  JKWrf'-2r*)  ( UV  fW-Zf*)* 


d* 


♦  2. 


.  JL 
<w 


zuc^-v*)^ 


\JA 


—  ot r 


(■)»*♦  r-1  ♦  o^-a^ry 

^UA^ir^V^  [  ■■■  1  „  __i _ J 

(V\ A  A  2yr)(VV fV*- ** f)’  V-lxf  )  (*•  '♦»V'-2*)J 


d» 


!l±[i-  -  fc*. 

C  ~  ~  A 


«  ^ 

-  t>,  vw' U')  e^J  cAim. 


'g  *  o 

^  **  £*+4)»,  ^  ^ 

*125:  j  T*(*0  *(».*)  ^Att^V  C,  (yi-M^e!^ 

*  „  T)( 

+  ^  f  |a.^> 

-  oi^v 


->7 


$r-$b 
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Tr- f*  .  *  ■»•  /  *  -*K  v>* 

S  ^  — r -  X  fwi-VT,  («►)«-  «!**»  -  **- 


-v^W 

X(*»r).  T,  c<.u* 


,  %[\fw»r).  7. (»••) -wW.  ,  T,^r)  T.(^*) 

-  >p*-  Jr — - —  «_  oIua  v  y»  — —— — *,>'c  *"** 

1  ..  r*r  1  yy'r 


1  1.0-0  _  3>lil  IM  f  * 

L  *l  »v.r  J  L 


-(x-w>W) 


_(cy)  [i-iy..  *■  iK'-4”1, 4  a^.>v-1')  ']  •**- 


with 


yv>  r 


i_  ’Lw1  -*■ 
1 


>  al  *. (  .,  .  _b — -  A  i^v.’^  <a  ^ 

^  i r  u'  L  (uWf'  J 

1  «• 

W1-  e*V  r7  -  2-^rtvi  »}> 


-  v* 
“  r  IT 


T  l 


Cr=>  *) 


£r-°) 

(r*°) 


-k|V)-w-)L*a*- »■"'*> 


A.  ].  35 


Trx*  Jo*  J, (W)  T, (mZ) ►,Wv'"Vw. 


cAaaa. 


witn 

«  .a 

cU^:  0  C*-*0) 

<0 

JTT  ^ 

- r,  •  (A. 4 

(kW)5*  T  T  1  ; 

|4 

,  +  "kL.y*  M~'): J-.M  I  V*A  *,***.  c,  (,-^,-K,ko 

J+mL 

—  T>,  ( 2- 4y.,  +  <- 

=  ?»(>"■). T,(».f).'3,[»)^-v^ 

E.  r  J.  *■  J. 

-  ^  V-J  ¥»rl^  l^C^l  (C,  t")  ^1^,-wi.):  tx‘a; 

v(A/X1'V'vlw“^  |T,i’)(i-S)~MyJ**A-VJ-k '  J)‘^“ 

with 

1  A- Cr-°) 


i!  V  *  K^!4 

*»  [fUxfv.wyM^'' 


1.  i,  ^ 


C^H 


-u-.  -  *£.>•> *  j.m.xmka  ufcJ;c,i«w.u^-T),(4-a)ia]jv^ 

E.  v*»  L  J 
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U-«  — - —  — -  «,  oUm. 

E. 


(-*+  hO 


J,  (»r).  T,[^0 


4i4i_S  tresses  _and_disglace(T«nts_in_the_gther_la^ers . 

The  relations  for  the  stresses  and  displacements  are  given  in  paragraph  4.1. 
The  terms  A^e"12,  B^m^e"1112,  (^me1112  and  D^me*"12  are  to  be  replaced 

In  the  second  layer  (H^h^Hj  +  H^)  by 


In  the  third  layer  (H^  +  h  4,  Hi  +  H2  +  H3) 

In  the  fourth  layer  (Hj  +  ^  h  ^  Hj  +  ^  +  +  H^)  by 


1*  * a 


,  -  ^V*»k  -  %') 

IS**  )•«- 

..  -,x  -c*’1'-*) 

(3>4*  >a 


A«“\ 
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ALGEBRAICAL  ANALYSIS  OF  ISOTROPIC  LAYERED  SYSTEMS  WITH  FIXED  BOTTOM 
AND  FULL  SLIP  CONDITIONS  AT  THE  TWO  FIRST  INTERFACES 

SYMBOLS 
We  write 

Ai  S  A;*.1  tW  Cic  C«~ 

c  _  E.  6**^ 

r,  »  - - 

k,  *  ^jC^KO 

L,  « 

E* 


H  m.  VW  H  , 

>y  «.  (.  H  ,  *  W  O 

X  *  W»  C  VI,  V  H,  *  »>) 

V  .  *.(».♦»*♦»>♦  Vu) 

where  Hj,  H^,  H^  are  the  thicknesses  of  the  successive  layers. 
The  index  1  applies  to  the  surface  layer. 

We  shall  successively  analyse  a  two  layered,  a  three  layered  and 
a  four  layered  structure. 


2A,«  4*  C\  (•«-/)*.)  - 

2»>,  •  A  +  C,  -  C-l-4^.) 

Adding  and  substracting  the  first  two  conditions  at  the  interface, 
we  obtain 

ZA,,.*-  C, 

lfcC’  -C,»*  *  *  l*)1"**  l*1] 
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where 


*  A|>*  f  -  C,[>» -2^ »-»)<•  4-  Di  (i-lyi**)  t* 

and,  taken  in  account  the  bottom  conditions, 

[ai]*  c***.1^  ')  &***  ♦  d.»* 


We  replace  Aj  and  by  their  values  in  function  of  Cj  and  Dj  and 
solve  the  system 

c**-  t^OLt1*-)****  *3  -  c *♦*»)«"*  »*3y 

J>^*e  «y 

2(4-2**)^-%^  .  „ 

We  combine  the  Xtt~  condition  with  the  w-condition  so  that 

tU"*  .  r, .  lc-'-hO  i>tf' 

and  with 


Fts  F,  L±±i- 
C«-h*) 


C,*-  ♦  .  Dit- 


We  replace  and  by  their  values  in  function  of  Br,  and  Do 


lA0I^+  e4>rJ  -  2(*+x)e*  ♦  2  (A-jt)l 

We  take  the  ta  -condition  and  solve  the  system 


-  3* 


V 


-  y 


B,r  [. ,  ;*-»]  „  ^  r  . »,) *  ju 


r  *  lAflLC**-)*1*  ~  O*2*)**  i*' 

*  — *  -  ■  -  — ■— 

V, 

l  -X  _  tAx]^***)^-  *■**]  »  (4-2>H'**-  «."* 

?, 

A,***  -lII/i]  D.*'} 

b,  «  2y>,  v  +  ^.(a-2^0  Cc.^V 

J),  *  /!_  a\A,*’V*  4-  C^^OCc-^K* 
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Chapter  2.  The  three  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  second  interface  (z  =  H-  +  h^): 

y  &jeY  _  C*  [a- «■  ♦  (*-  e 

***  4  -  C*f4.^§.y)«Y  * 

T-rj:  .  ^je'V  4  z(2)*i-y)  «-~7  =  o 

-  V‘y  +  CjC >H^)cy  4- 

>v  :  Ajt^-  -  c2  (.2- -  J)a.(.2--4|«.i*v)«-~ys 

K,  [A*«.y-  bj«TY  -  C>  C2-4|»>-v)tS/  - 

Boundary  conditions  at  the  bottom  (z  =  +  H,,  +  H^): 

*:  Aj**  -  bjO*  -  C3  -  *»( + 

C)S« 

2.  Resolution  of  the  system  of  10  boundary  equations. 

Adding  and  substracting  the  surface  conditions,  we  obtain 

2  A,  e  a  ♦  C, 

lib,  «  a  *  C,  - 

We  add  the  first  two  conditions  at  the  first  interface 

1byK  -  C,  4  $,«rxr  [A»] 

2b,e*  -  C.,t*  v  J>,  ■*■  2x)Ox  *  p,] 

where 


^  A2J  a  Az*X4  bjC*-  C?  4.  J>2  (>1-2|.»4a)C  * 
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We  replace  Aj  and  BIj  by  their  values  in  function  of  and  D,  and 
solve  the  system 

x  [a*]Tcz*-0C2%  a]  -  (*+i*)C*  ^Jl1 

C|»  a  -  -  -  —  — ■  -■  »■  ■  ■ 

2 ^-r2x*)efa-  e4*  - 

lAl3[C^2M)e%-e4y]  4-  (V-2Q£^-  e* 

2  (4 +2**)^  -  e4*  -  ^ 

We  combine  the  T^-conditions  with  the  w-condition  so  that 

ZC-l-^OQt1  ♦  2(4- |a.)  D,K*  a  F,  C'lft*  +  2(-1-}*i)  J)t£*j 

and  with 

r2  =  Tt  (A-y^ 

(«')*>) 

C,t*  4  ]>,£>.  ?x[Cztx  4  D***] 

We  replace  C1  and  D1  by  their  values  in  function  of  A2,  B2>  C2  and  D2 

z(>n-  2.x1-)  «.-**_  «T^V  _  * 

Together  with  the  1  ^-condition,  we  obtain  then  the  system 

-  &X4.  X  4-  Cx  (2)*!  +  Dx(_Z|**.*.x}«rX  =■  O 

Al*X  ♦  W*  -  Cl  4  *3>x  4  )x/*-h>t+  X  -Ss)CX»  _  *» 

.  fc, 

where 

ft-i  =.  7(x-x)  <  il!  '  a(44X)x'X 

V,  o  -  £4*  _  , 

a, 


We  soVe  the  system  for  A2  and  B2 

Z^z9  —  «•  +■  ~  2-x  +•  R'i)  e 

2?)jS  —  -?i  c*  +•  Ci  (  A4  —  X>2  (-'I- +  2-X  -  ^  s) 

f^i 

We  add  and  substract  the  first  two  conditions  at  the  second  interface 
2  Alt?  -  Ol(;>»-^1.2y)ey  *  2>2%y=  LA0 

2.|>xC^  -  Ci*.*  4-  V t. (_A- Ay-i.  f  2.y}  -  ^A“iJ 

where 

[A*]*  A*tv-  B>j4V-  cb  4-  [*-ly*+y)e  v 

and,  taken  in  account  the  bottom  conditions, 

[A>]  a  *  e^*  E>3.*y  *  ][  (x-Zfs+'f)  +-[Z-^5+l)«-  JPje7 

We  replace  A2  and  B2  by  their  values  in  function  of  C2  and  D2  and 
solve  the  system 

c2  y  +  **-0*^y  »(*-*»)]  »  *-**f>*  }] 

\.A2»3^(X4ly-2.x  v  (y|-2v  4-2x)^^-x;J 

We  combine  the  x  -conditions  with  the  w-conditions  so  that 
rz 

2  (A-yt.)  DZJ  *  z  ( 4 -  J>2%:7  a  V,  jl  P»tVJ 

and  wi th 
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We  replace  C2  and  D2  by  their  values  in  function  of  and 

(Zy~7*  •*•**) **■{'«-*%) *■  TrD44'/'^1^  7 


2^**  A  X*  x7_  -  [a-  R%) 

-i  Vct  t>>«tv 

Together  with  the  ^-condition ,  we  obtain  then  the  system 

^}iV  V^>|^  *.  [_ C* *  A)**  *  *•) *  *" 

>J  v  VV[(*'*H>  y)v  ij  *  - 


where 


..xNj-x)  .  ,  „  N  — 4(^  ->) 

§,=  a  [Zy- 2*  1-  *»H  v  fvt- R>»^  -  (-i^^3»)e. 

Gfci_=  ^  +  y-*)e*V“X>  -  £*-y*x)*J^  ^ 

V2  *  2  \  2*  x)-  f^+ f-i-R0 

*6-  ^1? 

We  solve  the  system 

-2.(a-y  ) 

«,  -V  4  C^>-y) 

9  m  _  ■  .  ■  —  ■■  --  1  '■ 


V’  =  *i. 


e  -  1 


7jsr  *  i*'  2y  V 


-2(*-y ; 


4.  ^ 
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3.  Relations  for  the  parameters^ 

We  have 

CiJ  s  — - ■ — - - - 

V* 

Da?.  MB” 

V* 

so  that 

W  -  I2-y  •f*-’0 

[^30£(4v2v_1x  yf »  [-1-2^ J 

Va 

We  have 

2 Aj*?  _  Ci  f-l-  -3-y)  v 

so  that 

/W»  i-^A!.]  *  -  p2.->. 

We  have 

2Tb* »  -  +.  Cz[*n-R>)t  -  Dz  A >«»  *-  2*  -  R >) 

“i 

so  that 

^Z.jC*  9  -£•  j-  —•  4  Cx*7  -  3)*e  •*•  2*  -  R*)J 

We  have 

LAi].{_C2x-i)c2-x  mJ  -  (^+  2*)«.  * 

J),e  s  - - - - - 

V, 

za,«t  -  c,  *  ®.**  =  l^i] 
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r 

» 


so  that 

A,**.  -1-  [LaO  +  C,**  C^-^v.-Zx) 

Using  the  surface  conditions,  we  have  finally 

=■  Zy*.,  ^  A,C  *  •  Ay,  2f*  .y  c 

3,*  >1  -  ZAt*.  **  V  C,«T  (>1-^0  *-X- 


A  2.11 


Chapter  3.  The  four  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  third  interface  (z  =  Hj  +  H2  +  H3): 

;  A>t.X  4  b}  t*  _  Dj  2^*. »-'*•)  •>*"  4- 

A*ex  v  W*  _  C*  (4-2 D4  f'l-'ZH  -*)  ^ 

Tjx :  A^«.x  -  bjc  *'  4  t  *■  ~ 

AAt^-bAe^  v  *•  D 

W  A*tx-  (_2 -  *■)  *•  —  (^2- +  **-) «-  m 

L  [Aa«*-  bAeX  -  Ca(2-Ah-»UX  t2-4hA**K*J 

U.  .  A^tx  V  b*4fc  v  c*  (/»♦*)** 

L  ^Aa«*  v  6a  t’1  v  Ga  (-4+--X.J  i*~  -  Da  e  ^ 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  +  H^): 

\v  r  o 

Aa^-  bA^  -  Ca  (2.-A^A-l-VV  ~  3>2(2-AH  W)th=„ 


Q4  -  o 
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2.  Resolution  of  the  system  of  16  equations. 

In  the  equations  of  the  conditions  at  the  third  interface,  is 
replaced  by  its  value  taken  from  the  fixed  bottom  condition: 

We  write  the  conditions  at  the  third  interface  in  matrixform 

B»)T.  M4  1 t>*  3>V)t 

We  invert  Mg 


where 


A  -X 

*  -  -t— - -•* 

A[a-  }m) 


I  i  {i-4)n+y) 

-w  »  (2-AIva*^ 

e,  1 

l  loah,*) 
\  L  L(A-^a*0 


-Y{Z-Ay*x) 


*t*-v-*) 


-^5*62 

*^5261  ^  -  C  ^ 


A2.  13 


vu,.vw«- 


r  * 

Q  O 
¥  ¥ 
i  «  O 


4-  ¥ 

\ 


^52  *^C2 


^52  6 » 


1-.  l-> 


Cj  o 


v  *(L.-'*0 


o 

I 

0 


o 


1  o  O 


L, »  -  a  (>~2}a%)  -  ( «*  ■»  L) 


Lx*  -  iL  2^-*)  —  1-} 

v».  i.^i.f-.-A^.)- u{»-«n)] 


e  H?2f2 


We  consider  now  the  boundary  conditions  at  the  surface. 
Adding  and  substracting  the  surface  conditions,  we  obtain 


A2.  14 


2 A, »  *  +• 

‘lfc,  =  >«v  c,  -  d. 

Adding  and  substracting  then  the  first  two  conditions  at  the  first 
interface,  we  obtain 

2  VC"  -  C,.*  ♦  t,  *'  *  1A0 


where 

[ai]=  A,%.x  +  K^x-C,^“^'x)e'X  +  M'1'2^**)** 

We  replace  Aj  and  by  their  values  in  function  of  and  and 
solve  the  system 

C  «.*  ^Ai^(3x-0  c2v  +  ~  1*)  C*  +  £.** 

3  e’a  ***3  *  O'*'2*)  -* 

-  -l 

We  combine  the  X  -conditions  with  the  w-condition  so  that 
rz 

lQl-y*,)  -V  3  (vt-^.)D,C,<  S  P  P?e  J 

and  with  T,  =  F 

c*~y  0 

C,  e.*  a-  3,  e X  =  F-i  v  J 

We  replace  Cj  and  Dj  by  their  values  in  function  of  A2,  R2,  C2  and  D2 

r  ^  ^  _*j 

If-If  ax^e2?  -  >» 
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Together  with  the  t  r2-  condition,  we  obtain  then  the  system 

A*4'**  b7tTY  -  C*  (>l-2^t~x*  v  D,^-2)»t-»x-R.>)]t"’Xs  - 

where 

ft,  *  a  +  Ax£x*_£a* 

aC^K^-af^Jcf 

V,»  2(a+1x’)£*.  Z**-  -t 

R*  *  ^  V< 

R, 

We  notice  that  for  m  =  «o 

Rx=  <3  V\-  -  *  R*a-Tv, 

We  solve  the  system  for  A2  and  B2 

»  -  ^-t.*  +  Ci  (-*+-&  a)*.'1'  _  ^a+-2*  -Ra^ 

Adding  and  substracting  now  the  first  two  conditions  at  the  second 
interface,  we  obtain 

2A,eV-  Cz  fs-Ayi-ly)**  ♦  Vy  =  \A&3 
+  Tl^  2yJ  ^  A  jJ 

where 

s  A$*.y  +  \b^y  _  C*  (-'i-2)*.*-y)e/  *  J,  y)e~y 


We  replace  A?  and  B2  by  their  values 
and  solve  the  system. 


in  function  of  C2  and  D2 


A2.  16 


c  y  lAflLfcy***  *  77  -  «-*(y 

«C^<| -  C" p*)^"-9-  H*0- 
^-T_  l^0L^v2y -2y+**)*^r ‘*}-  !]  V  ■p7  [•'^  -  (^-2yv2»)«.  y  J 

2^+2bH*>',('r'°  ^^(r)^',<-  oa^41'*”0-  i"'-84) 

We  combine  the  xrz~  conditions  with  the  w-condition  so  that 

!{»-(. t)  V,^=  w  [2(h-|-»)C»«>V  *  *C 

and  with 

W,s  wC^±2l 
C^-hO 

Cx^  v  h*-~y  =  *•  \.^y  +  2>>*V] 

We  replace  C2  sod  D2  by  their  values  in  function  of  A^»  B^»  and 

2^i{^-»)7^i7^K  X)-  (^"A0 

«  V,  ^C.*ey4- 

Together  with  the  Trz-condition,  we  obtain  then  the  system 

A,eV-  W*  C>(1-V*'V)?  - 

v  M^'V  =  ■  7oT 

where 

r  \  -^('f  -*} 

0,.  1  (1,-1.  ►  »,)eUv'T  l— *•*■)  -  R*> 
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$»  -  *i3- 

We  notice  that  for  m  =  * 
Q(  a  ^ 


We  write  the  system  in  matrix  form 


\  o 


)  o 


-[A-2y.b-y  vVb*) 


vy  -0^ 


^a-y 


yy.  m*,  *  *v 


mt]  i*>  *»  c»  *>f’  t-v,  ") 


We  replace  the  matrix  (A^  Bj  Cj  in  function  of  and 

-i  T 

►  N.  IT  "lY  -l^V-x)  rs*  4-  Mf2t2 

"  n?2?’  J 

-  \ 


We  notice  that  M4i *M5261  =  0 
M41,K5262  =  0 
M42*M5161  =  0 
M42*M5162  =  0 
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We  call 


M41,M5161  =  N1 


M41,M5162  =  N2 


M42’M5261  =  N3 


M42,M5262  =  N4 


and  write 


K  ''  v  *  Z^-*\lyWb  ♦  ey 


VAJ  L*a  Da)T 


(  %  •  j 


that  we  transform  into 


L« 


-2-VJt-y)  -*(>*-) 

„  +  e-  Nx  ->  e. 


rV^T 


-  °  ) 


^  D>4ey) 

T 


We  write 


1V„  i-  *-  IV:  t  (  V$  s 


a. ««  «n 

Ox.  pui 


rv* . 


\>n  w 

W.  ^>3 » 


where  all  the  terms  a.,  converge  when  m  =  o* 

* 


We  develop  the  matrix  equation 


(a.n  v  t y  +  (fta  v  ki».)  Da  t~y  <»  4  l*-)**,)  li¬ 

ft, 


(*a,  l»i«j  +  [»u  ♦  Wj  )  D*  e. 


-y 


e  o 


We  solve  the  system 


x 
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®<«v=  •<('<'('*)  ^ 


where 


R-tz  +•  i*z 

~?T 

Ali  4-  ill 

“ vT" 


V»,=  (  an  ►  lsi*)Cc'J-iV  J»40  'bii')(nj.i  +  in_) 

—  ^cxi,  4-  Vn)  cxzz  »  £c*ti  ♦  in}  ex*-'  4  b«- hiz  -  ■  ^2  i  +■  <*->•  ^2-1  -  °->i  •1)1> 

m  ^ ii * ^22  "  ^*12 *^21  con^a''r's 

and  has  to  be  devel loped  in  close  form 


The  term  b 1 1 * b22  "  b12'b21  contains  1 i near  functions  of  the  variables 


W  W  | 

1  ° 

i  o  4-2.|*j+y-&> 

f  0 

L, 

0 

L*  ^2(L,-U) 

1 

.  iii  Bn 

l-  ■ 

.1  o  1 

o 

l. 

o 

Li 

\>n»  1~| 

W  =  U^-atJ-.-U)  4-  W  [A-2^y,+  y -$0 

Vii  *  -l-i 

bi»s  -tj-zO-.-u)  4-  uC2-h>-y) 

Wfeii-W-W  =  x(i.~uz)  4-  L,(xh»-y) 

4^  valL,-U) 

=  L,U 

The  linear  functions  of  the  variables  have  disappeared  <•  ••  • 
the  numerators  of  B4<e‘y  and  D^.e’y  tend  bott  to  rf^  . 
factor  Q^,  and  the  denominator  tends  to  a  cors’ar* 

l'h^  <5,  L,  L,  s  (*♦*-'  -  - 


and  finally 


n  'y  -f*  ^ 
Bi  »  * 


% 

-  4  f 


© 


* 
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-  *1 

*i 

-  j: 

r  \ 

ft 
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3.  Values  of  the  parameters  A^,  D^. 

We  express  the  values  of  the  parameters  A^»  In  the  same  way  as 
explained  in  appendix  1. 

3.1.  Values  of  the  parameters  Ag,  Bg,  Cg,  Dg 

The  values  of  the  parameters  Ag,  Bg,  Cg  and  Dg  are  obtained  from  the 
relation 


(**  o.)1; , 


H5-202^j  (1*4  5-0T 


The  matrices  M526i  and  M5262  contain  nothing  but  zeros  in  their 
first  and  third  rows,  so  that  we  can  write 


(A,  O  c, 

1  L 


MflC2 

J 


(*>*  Mt 


ft  X  r.  4  \T  -d 

(A**  o  o  s  - - 


-O*) 
•  4. 

- 


4  *  tffisil  (  D/,ey) 

The  matrices  Mgl61  and  M5162  contain  nothing  but  zeros  in  their 
second  and  fourth  rows,  so  that  we  can  write 

(o  e>,  o  J,)T  »  -  -±— .  ]  £»H  M Wn.^  (K  MT 


•Y  JU?)r 
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3.2.  Values  of  the  parameters  A,,,  B?,  C2>  Dg. 

We  have 

W^t-  »  »■  . . .  — — 

Vi 

so  that 

Xc*  =. 

_  L''*3LC-’,1r1,vR')'^"-  (-'*I!>)^('',^J  '■ttL'1'  t'1-2''*'1’')4'  (7  J 

- 

where 

Ai,ty  4-  fej£y-  C»  (A-^-y)*^  +  (a-7y3>^)C^ 

=  (A,.1)-  v  .(C,t4> 

We  have  also  that 

2Ai«.v  -  Ct  *-  5  Lfti3 

so  that 

A**?*  |  ^A3]  4  Ci(A^-A^t-2y)t.y  -  \$iC*).iLy~^  j 

and  finally  we  have  that 

7>\>i  a  —  *•  Ci  [>+Rj)e*  -  "Dt  [s\-A  y«.x  v2x- 

n» 

so  that 

8,*’=  -f  1-^-  * (c,.?) 
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3.3.  Values  of  the  parameters  Aj  and  Cj. 

We  have  that 

C  iw  1>0-  +  *1  ~  v  c** 

V, 

xc* .  1*0  [(«♦*)  «-('■-»>  **"-  *  * 

v, 

where 

[A*]  *  AXS  *  C»  **  v  VM)  ** 

•  (Aiev)  e^"x^  4-  ftzey  -(ptcv).^-2^w)e^"*  5ze*. 

We  have  also  that 

iA,e<-  c,^-v-aOey  *  ^t',‘3 

so  that 

A,«.*^  LA0  v  ^  ~ 

The  values  of  Bj  and  Dj  are  obtained  from  the  surface  conditions 
A,  ■*  Q),  —  C|  (-4-2^.,^  4.  S,  (-4-  2^.,)  »  4 
A  ,  -  lb,  ht,  2^,  v  3>,  2|m,  =  « 

so  that 

b,«  2^,4  A,  A^»iC,  [a- 2j*,) 

«  *■  (A,*x)(.*a-^h') v  (C,*.v)  (.4-2^,) «. * 

3^  »  S\-  2  A, 

*  -1-  *2  (A,iy)  t"y  v  (C,  (>t-X^,)e“y 
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4.  Relations  for  the  stresses  and  the  displacements. 

The  relations  for  the  stresses  and  displacements  are  completely  the  same 
as  those  developped  in  appendix  1,  by  replacing  the  parameters 
by  their  adequate  values. 


Nevertheless,  there  is  a  problem  in  the  computation  of  the  vertical 
displacement:  its  value  at  the  origin  (m  =  0)  is  undeterminated. 

The  relation  for  the  vertical  deflection  at  the  surface  is  given  by 

which,  to  avoid  convlrgency  problems,  is  transformed  into 


W=.-^3W-  J» 


-  a^A.-f>*"2>A0c  J  ^ 


cLvm 


T*  (wf)  T,  (>*'  J 

*  l 

The  numerators  and  denominator  of  Aj  and  Cj  are  both  zero  for  m  =  0. 

To  eliminate  the  indetermination  we  should  develop  Aj  and  Cj  in  a 
Taylor  series.  Altough  this  is  theoretically  possible,  the  required 
computation  is  very  long  and  the  risks  of  introducing  errors,  in  doing 
so,  are  enormous.  Fortunately,  we  dispose  over  the  fixed  bottom  con¬ 
dition,  which,  for  m  =  0,  transforms  into 


Aa-Sa- 

The  w-conditions  at  the  other  interfaces  transform  into 

A-i  -  hr  ~C,(i-4ya)  -Da(2-4r)=  V  -»> 

A,-  lb,  -  c,  (2-4H  -*»  (2-A H-  f  (2-V) 


so  that,  for  m  =  0, 

A,  -  b,  -  C,  -D,  ro 

and  thus 

J_A,  =  *X 


so  that  the  problem  is  solved  without  any  difficulty. 


APPENDIX  3 
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Chapter  1.  The  two  layered  system. 

1.  The  boundary  conditions. 

Boundary  conditions  at  the  surface  (z  *  0): 

Tx:  A4*  +  +  Cxi.s>*  +  D,  ®  ' 

-  b,  ♦  s,  C\e J),  •  o 

Boundary  conditions  at  the  interface  (z  =  H^): 

.  A,  *  +  O,  +  +  Da 

■ _ 

W:  A,  -  4  **■ 

F*,  -  C*+^*)  bi  V  fa  (vni  |*a)  C zt  ^  J/1  fwa  f  J>aJ 

t*-  ••  (v>,+. y.\")  A ,  ■*•  Cv*«4V‘0  ^*e  4  ( 44  ^0  Ci  4  C^*h'^  ^,fc  — 

^  +  (»Xiy-i)‘bi  *- 4  C**4hOD»] 

Boundary  conditions  at  the  bottom  (z  =  H1  +  H2): 


*  •.  (*+  Jo.)  Ax  -  C^h1-)  &*«•  ^  "  +  «t  (vu.+  )*.l)  Cj.  -  fa  (u,  f  J«i  )  ^*^"2  o 
If  5a  >  >  Ci  *  o 

'b£{*~*\  ^Cjh^)  *  -b-*>Wrx) 

Xf  5i<-l  ^i«o 

C2-r,(y-,)^  +  ^  -2W»-0 

J*.[»la4)o) 
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2.  Resolution  of  the  system  of  6  boundary  equations. 

We  write  the  conditions  at  the  interface  in  matrix  form 


— 


T*  MT'ttx 

Ct>i  ®i3T 

*iOH 

-  r,(it.+H 

*> 

i»  (>‘44‘)*')  4- 

ii  (u  ■  ♦  |»>)  t 

-  * 
-*\t. 

-tt(n.4)*<) 

OH 

-  1 

,  \  x,r 
-r,(n,+^«)e 

ex** 

S,>-1 


I  , 

iw 


_  yx  4.  -(«♦*■)  CV"*) 

TuU*  OH 

L  C^vO 


If  >i 

,  a*Y% 

4  4  — *- —  «. 

Xv[n«4)>-0 

>n>t  -fiwOOH 

H-t  4-  - - -  t 

t*  ^vn.4^0 


'  4-  t 


-l/Ay-*) 


-A  t  *»e 


Fw[.)  4  ?*!-«  * 

L  Mwvif  J 
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*,x 

*i,e* 

R»e*  ' 

>1 

R>. 

2y,  (a-*.*) 

o  *.* 

ni4j  «. 

A,=  — ! -  \  RH-bi.+ 

^  a.  _ ! - bi  V  Ru  0l]  c-x 

2*,(a-v>i) 

^  ,  _ \ _  ^A*i-  &i-  *■  R»2 

2x,  (a-i»0 

X>(  a  '  <-  *  +  1 

2x,(a-v>^ 

We  write  the  surface  conditions  in  matrix  form 

Mo  LA,  B,  C,  3>Jt*  Or 

Ho.  MR  (bj.  ©Or*  *)T 


_y  -  fi* 

t,  t  e  > 


1  e 

-1 

*<«■ 

-  x. 

1 

Oln 

Out 

1 

Mo- 

HR  =. 

=  HA 

flu.  i 

i 

Q-n  = 

RuC*  t  Ru*.* 

ft- a  ’ 

-  %i*Ty 

■V  fUl*. 

*  *  »#t‘w  <- 

ftl\ 

r  Rut"* 

-  Ria«- 

*  v  f, 

XiRvi  *?' 

t>tv 

- 

-  Rix«.v  *■  i 

«  -X.x 
r, 

n  J.x 

-  X,  ('M7fc 
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We  develop  the  matrix  equations 

(X*  Iba  +  *«■  Di  “  0 

1  *.(*-»«*) - 

Oh  Alt-MW 


3)2 


-2.2,  (>»->»'') 


All _ 

£t„  Dltl-Al  i  A*1 


C\|(  <*lx  -  (Xitft  2.1  s 


*?*.*'* 


ft,  [-2Mi**cr,x  -  (*-*)  ~  l*«0  W  i*] 

4-Rz,  |^2ft|2C>e,,X+  (**•*)  ^*t1W  +  (-*-*0***] 

+  R*  [(-*-r.)g«eV£2,,‘  -  (*>»)*i2i2t'*-2t. 

4.JU.  C  V  2J,  *»**••  •f,>]^ 


-  t 


.****.  V 


UniV’*  -  gw  £**  +  x.  *n  C****''-  ■>.  *'»•“* 


^x*  —  2^« 


+  j“«  ’  -  j">  Rni  «• 


The  numerators  contain  only  negative  exponents. 

The  denominator  contains  negative  exponents  and  the  constant 

(>*-$. )  ^  Ri\.  Rax  -  g»2- 
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3.  Relations  for  the  parameters. 

The  values  of  the  parameters  and  Cj  are  deduced  from  the  matrix 
equations 

-  —  — —  [  fin  bt  ♦  ^*1  ^*3 

2r,  £*-*.)  ^  ~ 

C,a  — - -  *-  &>* 

2/, 


The  values  of  the  parameters  and  Dj  are  deduced  from  the  surface 
conditions 

«  j-(  ^  A,«. *  +■ 

C-+-  *0 


A 


2^*r*  -r  HCitM 


OM 


C  *-*•) 
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Chapter  2.  The  three  layered  system. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  second  interface  (z  =  +  H2): 

tz:  +  i>3  «■  », 

Xnt:  At-  -  x,  5, 

tv  t  (<it|4.)A{  -  ^  ^  *■  X'lYnt*')**)  C»  -  Jx  (Hi*  )m.)  = 

y^4-  fca  *•  £>«•  ^^  (<<+•)*»)  Pi} 

Boundary  conditions  at  the  bottom  (z  =  H1  +  H2  +  H3): 

W;  (4+]*>)  -  (-t+- J>(v»**-y.>>‘)  0*  -  ^  *»« 


■*»>  '  C^«  o 

A,;(*'y)a 

2.  Resolution  of  the  system  of  10  boundary  equations. 

We  write  the  conditions  at  the  second  interface  in  matrix  form 


(At  &t  Ct  .Pi}  a  M>  .  C 
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~  *1  (wtvj*t)  *t 

-  *7 

r  >  Cl 

ft  ft*)*'-)  «■ 

St^M. ¥  -  tiJ'1"** 

“**t7  A 

-JiO»v*k*) 

{*  *  )*')  -  ) 

,.  ^I-*-**) 

~Ji(Viv*kO«- 

H  *»>  • 


-l  (*-•/) 

K/U>H»0  -f^Jt)(2-v) 

\  ¥*, 

*  *  (u*%)  e 

4.  J*0»**O  „-(>-*-*  Ol*-v) 

~ '~’j 

( »♦-)**> 

y,J 

W-  (^y-.lj ,  4  **(.»»<+')  V'1*'011'10 

L  f*^)1 

f »  <  * 

4,+ub  _(*♦/»)  (*-v) 

)  +  - — '- —  *. 

i  f-  e- 

)*>) 

-*  *  y*^‘y) 

*¥■*> 

_ )  *-  — ' - «. 

+t  J 

VCvv  ih  lH*»-hO  ]_- 1  ♦  ft2  *  ^  y)  3 

WLfr**WP  v  r— * 

I  kwh*-^)  L'  +  V\ 

5ln 

6a 

i». 
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> 

■ 

K»  | 

4  - 

—  1  Rn  A ♦  R.J.K  * 
-*.)  L 

0  rs  -*(>-*) 

*>*  V_AX  4- 

R.w  Da] 

B>i  *  -  \  An  ^  v  ft 2a_  ♦* 

2^,^-h.)  L 

Rxa  Dx]c* 

C.  e-i 

—  \  R>*  Axe. 

.  Rn.'bi  > 

Aa*  Da] 

2*.l 

£.-».)  L 

5,^- 

- -  \Rai  Aie^  4-  Rax  fc* 

4-  Ra>  Cat 

4  Ry>  DiAe*** 

2t, 

(<•->*•)  L 

J 

( A ,  1b. 

C,  5,3 

» 

MR.  Mfc.  C  e>* 

»0T 

4x.«x  (>i-*i, 

V>n 

bn 

MR- M^-  )>L,t' 

bixx* 

b*i 

b*x 

Kv" 

v  »** 

BAXt. 

bn  « 

«  iV  (***) 

».  Rnvizi  «■  4 

RivAvi  *•  ^ 

bn  • 

R»t  tin. 

RixlXxx 

R»%  An 

Am  Avx 

bit  * 

Rxi  tu 

Rxx  Hn 

Rxv  A»i 

Rxy  dlvii 

y>lA  a 

RiiChx 

Poi  ®X2 

Rl%  djL 

Riy  (bvx 

bai  * 

R»i  (in 

R)x  6ai 

R»i  6a  1 

Aa*  ®a» 

iw  * 

Rax  &»x 

R»i  ®X2 

R»a  ton. 

R§a  Ivx 

Vai  1 

»  Aa. 

Rai  4xi 

Ra»  flat 

R**  €*» 

W 

a  Rill  CDll 

Rai  ton 

Ras  A*l 

Rkv  <Dvx 
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^1 

i 


We  write  the  surface  conditions  in  matrix  form 

M,  (A.  b,  t.  3>,  )T-  o 

The  value  of  MQ  is  given  in  §  1.2 

Cv«i  Ait 

li«(A  .  b,  C.  2>,)T. - i -  (>*  P») 

At, St  [*-*>)  An  au 


*».- 

\  “V 

t>il  «-  + 

V  * 

Oil  «.  V 

V»|  Z** 

V  All  e 

«*«.  * 

w*  «■ 

iiit*  *■ 

v  J.* 

4  AU4 

Pi  24  = 

- 

1  -** 
4-  S,  0> ,«. 

»  S>* 

-  S.bj, it 

t»tt  * 

w* 

-  w* 

.  -S,x 
4  X.  i)2 1 

-  -T.  tat  e6* 

We  develop 

the  matrix 

equation 

<*«  b 

^  +  Pi 

iv  Do* 

A  r,*\  (a-  + 

ll) 

011  i  bj  *  D3  ** 

O 

and  solve 

the  system 

• 

('1- 

Vi,)  (x*hv) 

Ail 

Ad-All  -  Att-A  u 

-At.jiO 

-Hi  )  (n-  >»\. 

^  01X1 

Ai-Au  -  An-  Alt 

7* 

^>1*  Pvll- 

Avi 
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We  develop  the  denominator  In  closed  form 

jt  s%x 

All .  Alt  -  ®-Xv  Aa» 

•  bnbn.^x* £s'*  +■  [<*-/•)  [_b»t  bn  -  b*  b*t^ 

♦  (■itftj  bn.-  bn  b/a^  %  ^  *■  t***)  J^bx*  bn  -  bn  buj  ^J,x 

♦  f4-J»)|b«  b*i.- b*i  b*».J  v  2  J,  ^  bxi  bii  -  b»i  b*xj  C *  ZT>* 

bj|  bn.  -  b»  bxx  =■ 

][0n  ®xx  -  fin  -  Ah 

*-  \>  -  *b,  **>  -  *»  <u*3*  l*~,'V,lvw0 

♦■  j^0n  Ovi  -  b<ti  A»v  -  Rii  RwJ  ^t.  ^ 

>  [0*0*1-  ab,®ixJi>uA,*  - 

*■  [0*  ft*x  -  0*.  OtiJfRuA*  -  All  *2*J  c^"x)  9?i(y-?t) 

♦  J_Qb.W  -0*  «*x]LAib«»A  -  «'*»**] 

We  write 

[0«  On  -  On  ^ ~**  *  J*  ^  «\»a 

14.  4m  -  o».  i*,,};*' ^V”1  (v”) .  «,,» 

[C>\«  0*t  -  Oil}  *  ^>4 

|_0u  0>x  -  Ob»  OuJ*-1  ^  ^ s-  C^ab 

[$i*  $iu  -  V$*  Ou  ^  s  ^24 

[^Ibi'&lz.  -  VBaaJe'^  oj5ij 
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k„  bn  -  b,  k..  - 


{ 


^  n  fl  n.  -  ■ 

Ao*  R>4  -  flu  AaV)  +  ^ifc^AaaRi* 

4*^ia  ^AhRh  _  An  Au^  +•  0|»a^Ri»Ri4 

.  «l'»-’V*w‘,°  <a*. 

».,W-  k»  bM  .  ,('.-’■>  J.W-). 

ojix  Ry*.  —  A  it  A  tL^  4>  Ojfb  j^flfc.  A  »>  - 

^A*  Aji*  -  flu  Aja^  *  f^xb  ^A»»  R>»  - 

*■  A1*]^**- A>*  -  All  A**]}  +  ®l>*  j^fl»5  RlV  - 

W„)  /.W-0  Qht 

b„,  ■  V  -  b.  b*.  J--)  %*«v->. 


-  Ri.  Ai*j) 
-flu  Rw} 

-  A.*  R^3  ^ 


1 


Rn 

Aox  R»*J 
R»*  Ajz.]jlj 


a  «- 


)  0|n-  ^Pu.  Rn-  -  *.fl**]  v  ®» **LN  R,!>  ■  *”  *  *0 


l 


■e  «, 


*0|U.  jjfl*,  R,4  -  Art.  R^3  *  *t«l*4* 
v  rtjii.  '[jUa.  a, I,  -  A**^  ♦-  *|** \%i,b  fl»A 


-  An  flrtjJ 

-  Ai*  Aa4^  ^ 


ifl  A 


»  t  v  i-  .lvy)  /*ly~x) 

Kt-  b$i  -  bia  z  *■  «- 

|  it >i  i_Rii.flw  -  v  n»»  [***»*-  **'  *a%] 

♦»J|4  J^Alt  A>4  -  *t»»  *  *>**  [*»%'  *»x  **0 

4-*)Ui  [A**  R»4  -  Am *  ^IRsvAm-  R*%  A*0  |j 
(w) 

«  e  e 


V&R4 
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\  v  v  \ 

Ma  -  *u  •  «■  * 

|  UtJ  + 

*4*  -  R*I  RuJ  *  Hi*b-  8*2  Rjj} 

*  ^i*  ^*3-^v  -  a**.  **43  *  v^^*****  -  *>^**43  ^ 

>*  i»t-  >„>*»  => 

I  ®)iv  ^Ajw  Rm  -  R.I  5*1^  »■  °li> |_*SJ|  *■»>  -  A»»  *a>3 

*■^1*1^1  R>*-  R*»RvvJ  +■  A^xR'dJ 

♦  -  R>i  Xym]]  *•  ®)>y  R*v  -  Mi  ^ 

a  ^V*V'L^  ®fcc. 

V-*  **  J' *  eCv  "x)  e  x  ^  ~*^ 

[^OA.C***'*  V  ^-J.)&&i£fc£2*,x  *-  I'l+S’)**)!* 
4  (>•+  *•»)  ®  A  4  «T2j’  X  ♦  {+-J.)  l&A*  4  2LX,  OKq  e  *  e  ^  " 

Between  the  brackets,  the  denominator  has  next  constant  term 

O-JOQR*-  *  O'0  [Qz,<Ui  -  viA.^»][Kii  u*  -  W  JU <,] 

Va.  t*  ej,x  c(>’x)  y' 
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We  write  then 


v  Cl-*')  \\  • 

b; :  .  1  «•  *  >J 


»« 

-tv*)- 

A  A  -(v 

On  a 

nx  e.  e  4- 

All®*,!*  4 

*«,$>»«-  «•  ♦ 

L* 

R*  ^»1 

R»>$H 

ky&rt 

!•- 

:  Rl» 

ItX  lu 

R  2*  tt  i» 

Rix  Am 

*«■ 

:  Ai\  V» 

Au.  i|zz 

Rz*  Ojr 

Rix  Uvi 

a.  R>  C)ti 

R>z  $!> 

R JM  ^Vl 

V* 

2  Ay  ftrt 

A»  Oti 

K»i  On 

Raw  tvi 

VAl 

x  Kin 

Kyt  Wbn 

Am* 

Rxx  &M 

V'*i 

i  %ln 

Kyi  llu 

Avj  t»i 

Ryy  fivl 

A-  y  c  %. 


x  *,x  {y-x)  '  . 

e. 


Ovz,  a 

|*  -1*-J,x 
®x«-  e  - 

■»1  -jr,» 

•zi  «■ 

w'xz  S 

O  -V-A* 

bn«-  «• 

-A* 

Aue 

+  x, 

4  C>i-V,0  c^-^o 


\»  .y  3s,t 

V  X,  *»»  e 

»i  -r -2j»  > 

♦X,  **!.*  * 

*11 


I  -X 


-  x  (x-x)  • 

=  A  g,xx  [*-*.)  h-*0  *  e - - - 

f  fto-*)  V' 

=  As,t\  (4-v»i)C^-hi)  iii. 

V" 

Dj,;  *±L_ 

V' 

The  numerators  contain  only  negative  exponents;  the  denominator 
contains  a  constant  term  and  negative  exponents. 
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3.  Relations  for  the  parameters. 

The  parameters  Ag  ancl  Cg  can  immediately  be  deduced  from  Bg  and  Dg 

Ax«, - \ 

C'i  a  - 1 - 

2  v*a)  J 

The  same  relations  for  Bg  and  Dg  contain  positive  exponents  which 
have  to  be  eliminated. 


— ! -  V  Ux.  *■  fiir  D*1  e 

2Jv(^-Kv) 


\ 

At,t i 

■V,0  (vj-x> 

2xv(*-H‘) 

r 

-S,X  \J  -** 
«.  -  oMe 

»  )  -x  JU,t  .»  -X-. 

♦  x,  b»r*.  *  _  r.  l»xa«.  J 

zz 

-J.»  l’  -3>* 

t  -  *%\  «■ 

V'  — X  — 2j»x  i/  _*-i7 

vy.  t  JJ 

Dx*  —  ‘  • 

r 

|  *■  ' 

V*  "  Bn*)!*]*** 

♦  «iw  v  v 

♦  V*  -  ***  CI*A]  ^ 

The  factors  R^g  with  the  negative  exponent  -Sg(y  -  x),  which  does 
not  eliminate  against  the  positive  exponent  (y  -  x),  have  disappeared. 
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We  obtain  in  the  same  way 


2s, 

V'  ’ 

^  °|ia  +  +  ^n>  «i*a]  -e  *  tX,t 

*  £\  \k  ♦  ♦  Rta  °>»aJ  e  *' 

-S,  [RvOp**  *  &3i  *fa*  *■  ****>3*^  tT*  ±7J'* 

W.  JU»  ^,i,  4“  #4a  ®»1.4  4  Raj  *t»4]  t*  l| 


The  parameters  and  are  obtained  from  the  matrix  relations 
V  _  I  *..  A,  c^-50  +  R,t  b2  ♦  WM'r*)CS  *  R,a  Pi] 

2y,(A-h.)  L 


=  — 1 -  \  R^iAie^  *  R*i  hi  +  R»>t 

2r,  (*->..)  L 


s,tv*x) 


C4  * 


R>a  Dt^ 


The  values  of  Bj  and  D1  are  obtained  from  the  surface  conditions 
2>,  «  ~ —  \  r,  -  (a*  /i )  A.C*  -  2s,  C,  ex,xl 

Ji-»  1 

*1  “* 
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Chapter  3,  The  four  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  third  interface  (z  =  +  Hg  +  H3): 

Aj  f  c,  + 

% rx  *.  A>  -  v  -  x*  D*  "^8 

Aa  Z^~*^  -  1b  ,4  +■  J*  ^  -  S'a  D* 

w  *.  (<i4-)Ab)  A>  -  D*<t  ^  J 

I.w  [(4+)^)  Aa Z ( 44-^)  6a  v  Ta  (yt**-)**)  C\ e/4 h,  (vtA  )  Da J 

U_  :  (b*  *■  )«•»)  A*  V  \-  C*  +■  ^ 

Lu.  A**^”^*.  (bA*W»)  &*  *  (^*^a) 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  H2  +  H3  +  H^) : 


We  O 


^  +  y^A)  A*  -O+f**)  ^  v  Sa  (h*  ih*)^  ~f* 


CO 


ij  U>»  C^.-o 

C>»+H 


\\  U  <» 


Aa  »  0 

:*(*-*>  j-«H)  ^ 


JaWa^H 


2.  Expression  of  the  boundary  conditions  in  matrixform. 
2.1.  At  the  third  interface. 


In  the  equations  at  the  third  Interface,  A^,  or  C^,  is  replaced  by  its 
value  obtained  from  the  fixed  bottom  condition. 

Me  write  the  conditions  at  the  third  interface  in  matrixform 

We  invert  Mg 

t/Ufe,C,DOT>  Ms-  MC  O*  SU)7 


til- 


-5} 


If  S4>  1 


^  ^  -s»J*  y)  - 


i 


k- 


-2(V-*> 

1  +«- 

-20*1 

-it«- 

.  .. 


.  r  jaIwa+W*')1  :  l/,vr't 


u,pTr^;p.  —  j  u,y*v»»*4r*)  L"^  J 

1  U  C^*>  J 


s4<‘ 


-2X4  (H*) 

^  \  e. 

-2/4(KjO 

-  Jj,  4-  e  .5a 
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We  write  Ag,  B3,  C3  and  D3  In  function  of  and  D^. 


A  e  00^4 

•>(*-**) 

v  2  Mfta).  (*->) 

vy  -  — . - — . .  '  '—  . —  1  ■  a 

c  _  .2.  moo.  MtM.  Em  +-  «g 

2/*^/t-h  >) 

^  _  ^.Mff^Q.y)cQ»)  B>A  »  Hy^O-^cOOP4^|a>v) 


where  M5(1,j)  are  the  constants  In  Mg1. 
We  write 

1>j»  =  Z  Mt6a) 

so  that 


A*  \ 

IT,, 

?*  ' 

*. 

>1 

TUP‘y) 

c> 

**(*-*0 

^31 

\r 

,  h  1 

TW^ 

T«e’,lM)  J 

2.2.  At  the  second  Interface. 

We  write  the  conditions  at  the  second  Interface  In  matrlxform. 

Ms  (. Aj  fcx  Cx  m  (As  B,  Cj  Dj)T 

We  Invert  M, 

(Atfci  V)y.  Ma  (A^bjC*  &0T 
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n».  * 


Si(4+y*)  -  *1  (*t*>*)  s, 

*»e«rt*t’H')  »,(«.*►*)  «.t»-')  -.,  .M 

-St  - 1 


-Jv 

Si 

a**1^ 


-(«•*) 

I  *.b*h)iM 


I 

-I 

-W,f^h3) 


k,  *>  (***>*)*■  -V,  r>  (y>-^4  |m) 

k,  y)  ki  i ) 


We  write  A^.  B£,  C£  and  Dg  in  function  of  A^,  B^,  and  D3> 


_  S.Y'lsfcQ.  +■  -S  ^3 

+  a.H,f.,o.vu/*K*,(v'^»  »  .s 

2Si 

Bi-  VUQ'Qe^A*  +  £  rt*[2Q.m(iQ-g>3  ^Cy-y) 

2*1  f»-  »-v) 

+  v  £  tt»(2Q- Mafia)  D3  £(v-y) 

2*»  (a-  >»o 


02  .  ^  tt»f»i).  TIaO*)*  1-  ^n3 fo)  MaQi)  B* 

(a«v»0 


•£  W>(*\  v  •£  w*  (*0 

JU»(a-vu> 


2/?  (a-*?) 
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where  M3(i,j)  and  M4(i,j)  are  the  constants  in  and  M4. 
We  write 


*w»*) 


so  that 


^  a So  ^  ®*  1 

C,  *  dWM  <S„  0»C’,M  *»» 

l,  I  1  QMiM£>W-«) 

2.3.  At  the  fifst  interface. 

We  write  the  conditions  at  the  first  interface  in  matrixform. 

n*(/,tb.C.  D,)t=  V\2(Azft,C*I)0T 

We  invert  M1 

(.A.Tb.C. -D.)T^  CtT>0T 


a-’  X 


-S,  5'  - r' 

\  x 

e. 

5,  -S\L 

-r.  (vi^O 

(A+Y'}  *'  rv 

#  x  J.X 

•l"  *0I',+IM» 

.(44).,H5,X  «-**  S,e  ’ 

-MH 

-nly-x) 


^  -fi  ^  (»»*♦)* ■*) 

M»»*+r)  ft  [^Y^n[y  v) 


We  write  A^,  B^,  and  in  function  of  A^,  B3,  C3  and  D3 


A,  * 


* 


tii^i)  Mips)  «rf|(y~yV2  *  xntb‘0-  M  jk 

2  X,  £»->»,) 

- -  -  .  -  ■  —  —  .  e, 

1J",  (*->\ «) 

T  ■■  —  -  ■  —  ■  ■  ■  i  - — ■  ■  *  £ 

2s,  1*4  ~)r, ,) 

25*i  ^ 

TS  51  M«(4'0  .T1»Cm)4L  ^“^2  +  ^  (4>  >  M*0O  G>7  ^j-rX 

i)f  __.  - - - ...  -  ■  ■  . « 

2X  | 

21  M«  (Ai ) .  M-,6‘-j)t ^^2  +  {A>).rh  04)  £>z  s  x 

■¥  m  ... - - 1  —I—.- . . . . . . . - .  . . .  - . -.i— .....  ■  - .  -e 

2.X,  (>*-*.■) 

where  M^(i ,j )  and  M2(i,j)  are  the  constants  in  M^1  and  Mg. 

We  write 


Rj*=  M.CjO-^0*) 


so  that 


A.  | 

ft, 

A 

C, 

*J\(*-*0 

*1 

u 

fc^V 

%u,S 

Rvi^* 

^>2 

W*?f 

W* 

Ax 

K 

Ci 

Di 


amaa.*  tou^ 

Wr1te  the  "ndU,“s  •*  the  surface  matr,xfe™ 

|  «*  -«  «*’  a 

L-  .«  (Ai6-c's°T- 

"*  have  then  follow^  systm  of  matr1x  ^ 

hr .  (A.ib.c,  d,)t=  of 

(A,fi>,C,  D,-)Ta  -  MR  /■  AiB»Ci  D2*)t 

[Azft,C2Dt)%  _J -  IiG|  (A*B3CM>3)T 

Jj* 

CA^CjD*)7*  - MP  [  1U)T 

i  **{<»-•»>) 

so  that 
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Wv(a-r,)C-r.Hr-M  '  t*  MT=  f*  o)T 

HI.  MR  =  r  **  "'*1 


1  OL1) 

CLJ.Z 

I  W 

W 

MQ.hPa  *>»• 

bn 

b*» 

biz 

i  b  4» 

b*z 

MI.hR.MQ.HPa 

r 

•  ** 

/»3.  25 


C21 


Da  =  -  Sx»  rar*  [*-*“0 


C„.  Cat  -  Ci*.  Oxi 
ca» 


c.«.c>i.  c.»2-  c*i 


3.1.  Determination  of  0^.022  "  C21 *c12 

The  expression  of  the  denominator  c^.C22  •  C21  *ci2  must  be  estab^shed 
in  complete  closeform  a 1  tough  regarding  the  exponentials. 


c„  as 

CLn-bu  4-  G-\z-  bai  ♦ 

Ota-  b}i 

r  ftu  ■  b«i 

C>x  « 

o«.  bva  4-  oui-  baa  + 

v  oim  .  b^a 

tat  * 

«H-  bu  4-  <*•!*»  *>2,  + 

Q»vKi 

V  OLlA  •  k*  ( 

Cia  s 

nil-  bn  *■  oiibu.  4 

.  cn.aa.b52 

4-  A.2A  •  b*Z 

c„-  cn  -  «,x-Ca<  « 


Cut  "  b>»a-b*') 

+-  (on,-<Xi-b-  a» 3.t*n)-(  Wbn  - 
f  (*n-b*A-  ciiA  w*0  •  *****  “  bn*bA») 

*  (*,*.  ai3  -  A>V***)-  (bat-bja  >  Wb,0 
V  (on.om  -  a, a .e^OO^W-  *»* 
f  ^**.*14-  m*.  «■»*)•  m^kl~ 


r  - «* 


CL 

0.12  * 


*.  V 


_  -r.x 
«»»•*  * 

_j»  x 


OlA  - 


L« 


tue 


►  R«t’”3 

v  1U*  «-f'xJ 

►  R*it 1 *J 


A 
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<v*. L*-. *"  -  *— *  * *•**■ <J 

■ «»-"  «. w"  *"] 

|V*  -V  *  A  *«£"'. 

[W  -W 


CiJV**' 


-(v-*-)  r.x  x 

Ct»|-Oi*X  -  cui-  *  *"  •♦'•*'• 

*«.  [-2^^'^’'  +  -  (*♦•)  wixJ 

* ju,  ^a*ner.£M  V  *  (>*-',)  R**] 

+  ft*  ^-rORn£a,£1/”f  -  (^*.)8«£2J,X  -2x>^fc,‘'t'  ’  ] 

_  -X  -/.X  1 

f  ft*,  +  t  s»-0  ***  v  2r,*wt  J 

-lv-*>  V  «/ 

o.v  tx-**  -  ®-iv  «  e-  •  * 

ft,,  ^-aft**^***'*  *(*»-•)  **»*  *  '  -  (s&'W***  ] 

»*»,  vC»*«l  8«£’"  *-t^**>3 

♦  *»  »•«*-* ’ J 

*tU,  ♦(*-•>»“  *  **WV“] 

Ji* 

an-  tM*  -  CM<4- <m.\  c  c  „iK-» 

„  t„-,  j 

„U.  *****"••■  J 
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Ru  ^-aRaa  C*«,*'v  f  _  (■*+*%)  R4»«“^ 

*  Raa  j^  a  R  »  *.“*  C*1*  *  (a+Ji)  RmaT **'*  *■  f-^-xOR^j  J 

'l** ^4 h**± **J 

*  Riii  j  C^rORo  a"2*  +■  (r*-')fci%  *  2f,  Rjj«rV','x3 


X,*  X 

®-ta*  ^24  -  ^1J|.  ft22  *  <-  ■  •»  • 


#,*  1-xKxt  CV1’’  v  (,'.-■>  (*».)**  iU] 

4  Six  [  ifu  «rV*  4  {»-»•)  *M  •T2”’  4  {*-'■>  **  J 

4-  R 92  ^  O-j,)  *N  e5**2*’-  (44 a)Rx4  t2*”  -  ir<  R**  * ’* **] 

+  f?AJ  J  (•HX,)*>*e2’'  +  U>-0  *M  4  J 

-Jii^-0  j,*  * 

0.1$.  a.  14  -  *»4-  **>»  «-  •’•*'• 

R,„  £-iti4eY*/'*  4  ( r, -.)«>*  i**;”'*  -  t44»')''« 

4  ft  i>  [  4  tvk  *  ’  C,rf  4  (^fl)  *”  t'2J’"  *  ^  *"  J  , 

4 1  „  [(,-„)  D« iU^'- «4'^ *'*”  -  2/'  %*  \'P 

[Ow,)«*e-J'  4  4UW'.*','J 

-1^-^)  Xi*  *  *  „ 

out-  -  »iv&n  »  e.  •«••*•• 

»„.  „„  -  „„.  m,  -  .-***£*>-»  .*V.  A,* 

»XA  -  «x,  -  A>* 

-x,(7-x)  x>*  x  l_. 

WV1.  Ol>  -  *.»V  » U*  *•  •*•*■'** 

ftrx.  ®ii4  -  ®(iA.  ftn  *  *•'  A24 


*1V  *24  -  *•*•  **% 


3*—l>  ■*»./•*.**.  A,i 
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v  $n- 
v  On-P* 

4-  OxvPxi 
+  ?ai 
V 

♦ 

4  $4l?ia 


*  ®,v?»aM"y>. 

4-  $rt^2 

+  *hi-hz 

4 

4  ($».  ?J1- 

♦  W^' 

4-  $*»•?» 


\  q,».^,e”tI'V) 

4  ($14  2 

■»  ^.Pai  J 

♦  ^J4.p4l3  e^  ^ 

*  Qm-P”41 
+  $*A-P« 


.  vvil"y). 

5  ^11-  ^ 

X  *  ]$2V  P'1 
>\ «  $*VP* 

*1  =  $*vP»1 

h\  = 

Ai  = 

-  ,  -l.-v)  -»(-> 

♦  !«,,«,*  -  *■•»){*>■'"  - .».(-) 
♦  Ttyv^*1  *  _  -Ix-y') 

» l*. t» * \  • 


\.?n.T»-^.wO«l<'y,i',“'Y) 

ivK“V  ;3a-^>^ 

>.1  »  l^.’n-  ?’0 


)P<»  * 
>'4 
>»* 
v» 

>** 
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W  V**-  *  eC**v) 

*^0\i0w  -  On  onj*  >*»  *  -  o* 

WW-WW. 

j  ]%*ki-t\iA*'y>n- +  o.a.oa.]-v»* 

+  [0irtfc>-  0»  #*i}-^*»  ♦  [$»>$*>  “  On  O**}  >*’ 

Ou.®**0' •>»* * ®»i'Ca*3  ^42J 

,  ls-0  (**v)  _M*~v) 

ta».  V*l  -  Wl.  *  €.  •  ® 

I  l^  Ov- *»•*»>]■ >»  *-t*1.0**- «»*>']>« 

♦i**.** .  «»..*»o->>*  *  i*-***-  «>“»’*>*  i 
.  .  .  .  Jtr* 

W  W  -  W-  ©*»  «•  «■ 

o»vO**3*v« v  ~  ^'  ®AS3  ^ 

4-  [Vv^4A  '  <&»■**»}  >J*  W  [^  ^lI  ***  )• 

>C  }* 

*|V-  ®A>-  ^0  >13  ^  ^V^i"  ^ 

^^A»A23 
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WW-Vn  k.,.  '>.(>■*) 

u  v«-  ta.fc*.  . 

W.W-  V. ■*  • 

WW  W.V„.  etv-’O.C-U».C.-v)  fa, 

Vn.  \>A2-  W- Wh.  C  '  •«.  •  «.  ■«- 

j 

Hi.  b42  -  >J1 *  *• 


The  value  of  the  denominator  is  then  finally 

y  r, y  N-*>  »*Y-*>  C*^> 

c„.  C«-  c„.c2,  «■  t  e  •  *  -  e.  •* 


[Au.  ^»u 


-tv- *) 


.  e. 


2^-vl  ^a>>(v-y) 


*  Ay&.  Sn*  e 

»  Am.  €r,(v’rt 
+  Al>.^e»^-') 

4  A  24  fc**  N  _ 

4  A*4.  ***  *  e  J 


The  term  A24‘B24  contains  the  constant 


$44  *  (*»“0  $4Z • 

|Qi4<ha-  QUi  -  ^3?20 


We  write 


c«- czz  -  £.»!•  c'ai® 


X  Jy-^.  e^'y>  C-,2 
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3.2.  Defcrminatlon  of  the  parameters  B4  and  C^. 


We  write 

a*,.  -  Hi.  «w  ♦  *.  *»  -  *•*•**] 

*V-  -w"  + 

.  U„. 

L  .  /\  0 . . 

Vn.  = 

a.  j_tyivSt 

5  j^a- 
i'u  r  9» 

Vi,,  -  [W 

}\t  i  ^(W  S* 


4 

♦■  $aa*  S> 

4  ^22- 

v  tya*Vv2 

f  tj*2  ^i 

*■  $H-  Si 

4- 

l  Oai-P** 

■*•  di*>-  Vy 

4-  QAl  .  Si 

*  ($4*.  JW 

♦  (Da*?*  ] 

f  $24  Sr*  J 

+  %A^>  J 

4  $2*4  • 

4  ($44  •  P*Z  J 


We  have  that 


$/,«  if,h  rilA-h^^-Vir)  (4-«»V  ■ 


c»a 


Cm.  c-xi  -  Cia-  £-*» 

c  Zl _ 

C„. Cxz-  Ca.  cat 
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C-ia  *  c*-a».  fen*  +  ai41^2 


„x  r,x  f*-v)  T  <Y~*)  '  \J 

*  «.’.  «" 

^-’! ft>,,.\,W  *  eJ,1^W-l>’» 

v  e-(v-^  y>5  4  !>«  J 

cai.  cx,.U  V  a-u-bi,  4.  «n  h,  4-  a».w 

.  fv-*)  >  «)  -Mv-*)  t  u»  n(y-^)  ,  w  n 

♦  *.'r  •o'aa.bx,  +  t  ©'iv*,>>  *  ^  b*1  I 

=  t*/.»  ,Cv-*>  t(*-W  .'.(cO 

r.2(y-,).„(y.,)  ^  Wii 

X  cN-»> ,>.tl->)„,vki#  „  e-Cv-y>  ^ ] 

tn  e».»  ='«■ 

»  I.<  N-xl  (*-V)  «'* t?'Vl  Cxi 

Cx,  S.  e  -  e  •  ft.  .1  •  *- 

>  c’** 

£>4«  %*,*>.**(*-*•)( :'*-*-0[->'-v*‘0  £l“~ 

cTT” 


The  numerators  contain  only  negative  exponents. 

The  denominator  contains  a  constant  and  negative  exponents. 


r 


i 
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4.  Values  of  the  parameters  A^,  D^. 


4.1.  Values  of  the  parameters  Ag,  Bg,  Cg,  Dg. 


The  values  of  the  parameters  are  obtained  from  the  matrix  equation  in  §  2.1. 


*»  1 

r* 

^7 

J 

4 

ft, 

ft,..1*-” 

Cj 

,  *»J 

ft* 

One  obtains  immediately  the  values  of  the  parameters  Ag  and  Cg 

/U.  — -  I  Pu-h*  *  ^*2  ^-**1 

— - —  [?>,**  * 


The  determination  of  the  values  of  the  parameters  Bg  and  Dg  need 
some  more  computation  to  insu?£  convergency. 


lC’-Y) 


x  "Pii-DaJ 


_ . _  |  P-4».  ■*.  P42.  D/,  I 

-^1  J 


Those  relations  contain  positive  exponents  which  much  disappear  to 
avoid  overflow  problems. 


Bi 


eC*-y). 

2  .  £,2 


^'v) .  8  r ,  Va  -  ?A2.  Da 


?U*  •  C\a 


‘] 


r 


[  i' 
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r 

F 


?*,.*>*  -  • 

iCN-'\-’'>('»-’)CTVvt'.i  *  *rt’r’’'  »«• b!l1] 

♦  e-b-’V*  •  Vai  ]  = 


0V5Y  Til 
+  e-l'l-’)  A*  [  W  to'«  - 

-M?„  Pn  v  ‘  **‘*0 

Mh,  [a,.  Pi. 


e 

-e 


4>u  T*'] 


=  .  ?%>*). 

eP  Pj,  p».  vn«-  .««  R*0 

T«l5».P»  **B  * 

.  A™)  ?fc2* 


x 


M.  35 


Po.>.  = 

^  4-  (D*v  ^2t2^a  V)+  Hvt.tito] 

-  «Tla‘V^  ?ax  ?n  +  ^'] 

s  £^y\  ?&*> 

'Pl\  'b\'L  -  Wa.  fcii  « 

e~(*r)  .  ?.,  ^Qw .T,, -W-'*!1-^  »  lj)4  j .  V,.e  2  J'^'V  +  <W  ?A.] 

-^.Pa^A  *W  P».  4  ^  P*J 

«  <■'  <*‘V>  PftjA 

The  positive  exponent  e^z  ^  can  now  be  eliminated 


£>»«  4y,Xa(>l-V»0^-V'2) 


L 


-M -j..M  „.tv  ve.t>  *  .-^V.a  pv]  -^ 


+  «.  -e. 

and  the  numerator  contains  again  only  negative  exponents. 

W  ft*  -  P«.t>A  = 


Pa,  \ 


4  i-W-iP.-lr.W-Vv*.  V«  »  e.'tV’1 

Wn'^i 

J\l  wl  -TMv-xl  \  tl  . 


-PAtU^'V^V,,*',,  h-",! 

^fW-^V’^’Vvv  v„  +  tl-','v)  *'*■  b’Ai  ]  - 

JW-Khlr^,  ^v>*.  b\.  -Pa.-  &}  *  €S‘M»'»3>'  V>'«-P«  b,1'3 

^1-\11'MI>'„,^?4i.bl41-  Pa.  -V».]  +£*,JyVi»^'>«-?«->‘‘'  3 
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^  vx  -  ^A1‘  b  "  * 

.  fe.Pu.»,vt»C|M^M] 

_e-M»-v1 .  ®„.Px,  *  vll->  J 


=  P6*i 


P^I^XX  -  ^Ax.  ^X»  = 

*  f>4l  «». ?». 

_e-j.(*-v)  *■  <5«^ 


♦  o^.Pw^V^J 
+  Pj’  J 


=  e-**t»-*>  TbA5 


?i!i.  b^x-  = 

pAl  i>,1£1,,'y) 


*  ^>1-  Pl2 
*-  1^32-  P2* 


+  $33 •  Tvx 

4-  $>*•  ?»» 


1 


s-Ma-v>  ^  .%x;2,a'y)  fUlAx-Pw-  *  *ie  e  ^ 


=•  e 


Yft44 


and  the  positive  exponent  e  3 


s,(z-y) 


can  be  eliminated 


Pj,*  4s.  *■*  • 

^tr-9-fch-^  PbA1  , 

,.  »-\W  ?W  * 


~ 1  L»  -ix 
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4.2.  Values  of  the  parameters  A2,  B^,  C2,  D^. 

The  values  of  the  parameters  A2  and  C2  are  immediately  obtained  from 
the  matrix  equation  in  §  2.2 


The  values  of  the  parameters  B2  and  D2  are  obtained  from  next 
matrix  equation  (§  3): 


^11 

Vn 

V>?i 

V>2 

bin 

The  necessity  of  converqency  needs  again  some  more  computation. 

. — i - -  -  *“-*0 

^S,  *•*(*->»)  (*-*0 


3  ?.(»-" 

>1 


1  \  *- 

-  C»-WJ 
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c'aa.  tat  -  c'n.  tax  « 

Vn  v  .v„ 

V  b'„  „  C^U  H  !>*• 

.  W.  v 

V  t,'‘l  ] 

„  kVWj 

+.-Htv-,)c'„  ^.t,*  -v„.w*3 

4  J_t>'M.\m-  b'v.^O 

\  b!*, .  luj 

-2h-y)  -rtf^-y)  \  Of-*)  tl 

at  -  «.  -  8.  «ave  ■  »tx 

4  eT^  .tjK  bx* 

4  e'^  o.'a4  .  "x^-  bz A 

B<4  4  C"‘W‘“l>’a4.bz>  - (rt/-’,)t'^-x)^,.'B^ 

d-\2 

The  positive  exponent,  e(y“x>  included  in  b21  and  b„,  has  disappeared. 


c'aa-t>it  -  ckv  ^AT.  = 

V*1’"’  [W*-  W.-V.] 

4  [V„.V*,  -  V».  WJ 

v  *14  [VA1.W  -  Vtl.W] 
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1>4\  -  " 

_  e.  . t  ©*3u  c.  • 

«>«  .^ly"')- 

.  e-W-^  e^W-’V*,  ■  es‘(r")  Bm 

%  x  .  )  -(>r*)-JiCy-*)  o  l 

T^,  .  '  £*W-«>  «'*  *»*  *  "  •*•  ' 

CT,x  L 

The  positive  exponent,  es2^y’x^  included  in  b^j  and  b42»  has  disappeared. 

The  numerators  of  and  02  contain  only  neoative  exponents  because  of 
the  presence  of  the  factors  a^,  a^.  a^  and  a^. 

4.3,  Values  of  the  parameters  A^,  Bj,  Cj,  Dj. 

The  values  of  the  parameters  A j  and  are  immediately  obtained  from 
the  matrix  equation  in  §  2.3. 


A,- 

C,  = 


j* _  T  v  •  C* 

It,  (a-*"') 


_j? _ B*  *•  8«t 


>  •  t-2 


■v 


The  values  of  the  parameters  and  are  obtained  from  the  boundary 
conditions  at  the  surface 


A  ,£*  -  2j\  C.e  J 

2A,ty-  f^+-J'0U|  *"  "I 


A4.  14 


2.1.  At  the  third  Interface, 

In  the  equations  at  the  third  Interface,  A4,  or  C4,  Is  replaced  by  Its 
value  obtained  from  the  fixed  bottom  condition. 

We  write  the  conditions  *t  the  third  Inter**#.*  4* _ _ 


APPENDIX  4 


ALGEBRAICAL  ANALYSIS  OF  ANISOTROPIC  LAYERED  SYSTEMS  WITH  FIXED  BOTTOM 
AND  FULL  SLIP  CONDITIONS  AT  THE  TWO  FIRST  INTERFACES 


SYMBOLS 
We  write 

At/-  -v  w  C  *.♦  —  ^0  A 

A; L  J  -  Ai 

t  *  -viaLH,  *Ri>  •  H;.*) 

(w  Mi  V  J  t  *  Bi 

C;»v»  njj;  (Vw  *)*•)*  -  c< 

3)i  W.'-M-.  S:  (*;  v  y*.', )  e’'*1'  (  H  ^  *•*  ^  . 


F,.  £jl 

£1 


V,  T 


L,,  £l 


X  »  KV»  H  , 

y  »  v*  (_V) ,  »■  H».) 

' Z  3  vn  C  W ,  V  VJj.  ♦-  W3  ^ 

V  A  Hi  +  H3  v  B*) 

where  Hp  H^,  Hj,  H4  are  the  thicknesses  of  the  successive  layers. 
The  index  1  applies  to  the  surface  layer. 

We  shall  successively  analyse  a  two  layered,  a  three  layered  and 
a  four  layered  system. 


A. 


A  4.2 


Chapter  1.  The  two  layered  system. 

1.  The  boundary  conditions. 

Boundary  conditions  at  the  surface  (z  =  0): 


Tx.y  A."  +  B,  +  4-3),  =  . 

1f£  *©  Ait*  -  b,  -r,3)|S« 


Boundary  conditions  at  the  interface  (z  =  H^): 
Tz:  A,*b,C%  C,  +  3 ♦ 


t-2e 


Dt 


trx :  A  ,  -  ♦  JiC,  -  Ji  i), t 

na«-  -  *  ftCt»  7  .  Di  r© 

w  :  A ,  -  «■  s,  Cw.+y..)  C,  -  S,  C^.+hO  3,ls‘\ 

Axt  ^  b»  4  Sv(v»i>)*t)  Cl*.  ^  Di 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  H^): 

Wa«  C'^^l*7')  Ai  -  |*i)  Cj  -  +  ^  *® 


Lf.  sxy )  Ctco 

Aw1'**'0.  *  i ±^±±Th^-''‘->'h-*) 

r»  C 

2.  Resolution  of  the  system  of  6  boundary  equations. 

We  add  and  substract  the  surface  conditions 

2A,c*  -  *  -  ^-j.)  J>, 

lb,  «  ^  -  ^4-X.)  ST,****  -  [A+i,)3>> 


A  4.3 


We  add  and  substract  the  first  two  conditions  at  the  first  interface 


ZA,  ♦  O*,)  C,  +  0*0  D,**'*  c  [Aa] 

2&,*x  V  )C,  ♦  [A  +  t, *  [/U] 


where,  taken  in  account  the  bottom  condi tion. 


LA3>  ■*"" 

if  <  ' 

lA1>  ‘  ]*»  -  t— 

We  replace  and  by  their  values  in  function  of  and 

v  (a-s.)  i.  J_e*'  _,*]  =  ^AiJ_  .* 

(A-M  C,  {_».»*  i"]  .  [-.«■)  3.  .f*]  *  [At]  - 

We  solve  the  system 


rtlmh 0  } 

r  . .  i  ’  J 


C,»  {  Oj,}3 


+  -  (***•) *s>  *  (^-*0*-  «•  3^ 


3>,=  |[A0[7,-t"X  -0*)*J,X  +  ( j-,)^*** * x** ^ 
+  ^2j-,  i*  e-r,<  -  O«r0*lx  ^O^]] 

v,.  .  oto,i«-1x*€2j'*]  v 


rvi  s  oa 


V,  •=•  ( -'-*•)* 


A  4.4 


We  transform  the  w-condition  utilizing  the  r  ^-conditions 

c,  -  *■»'*'*  ■  F'  W»*-»  0,] 


and  with 


a  '  t  — - - 

j,  0.-0 


cx  -  a,***  *  F*  i  )  -  *'l 

We  replace  C1  and  D1  by  their  values 


C,  -  »<e 

[  1*01 


<1-  «2x/X  /  \  —l* 

£/»-*,■)  -  [*  +  /.)£  v  (4+-Me  -  *  J 


.  .v,  .  -  -2*  -/.*  2r5'*  VI  - 

4-  \2*,t/  (/»-  t  )  *  2*  4  -  ^  J  J  ^ 


n  Vcv*r"(v~x)  -  »x] 


.  f  _  -Ji(lf-*)  ,1  Hi 

]Al]  -  *»  LCv*  .  l,j  =  -  — 


This  relation  becomes 


**>  » 

W**  ’)  1-*  -r^ r  ‘  *  ^D,*‘ 


ip  *i<< 

— - — :*  e  J 

u  n(viv^»o 

,  [>14  Pi«--*2. 


3°  ? 
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Together  with  the  T  -condition 

rz 


\[  » 

K  V  L  MwH .  ■>  ,t„ 

L  J  u*  H  J 

we  solve  the  system. 

\\  *»>  » 


*tCv*»V^O  -bf-*)  -Jvf'f-x) 

®x«  -  £i-  -T^To  ■* _ 1 _ 111 

*•  V2 


®x-  JSl  J! _ ll 

*■  v* 


H  n<> 

b,=  -  3x  - 

*  Va 

x  vyy  -I'f-x) 

Si  2i^hi *  _ ll 

*’  7t 

„  .  .  -2^7 -*) 


aO*H  -J'(y~*) 

-i  1  **  ^ 

(ni**') 


3.  Relations  for  the  parameters. 


c,  •  *  (*-'•)] 

♦  ^2x,C*  -  (h+J-.W1’*  V  (a-*.)  c1'*eJ'*3\  ”~— 

il  lAl]  -  (ava.)C,  - 

The  values  of  the  parameters  and  are  obtained  from  the 
surface  conditions 

b,-  -i -  1  X,  -  (4W.)  A,  a*  -  2r.O,  eJ'x] 

jr,-»  L 

D  -  _  J _  \  i  -  lA,e-y  -  (>+'•)  C,  C**  ] 


A  4.7 


Chapter  The  three  layered  structure. 
1.  Supplementary  boundary  conditions. 


Boundary  conditions  at  the  second  interface  (z  =  +  H^): 

**•  A»»  VCy'y,.Cl  *  >,  V  ^'y\  P. 

Xrx:  A*  -  bze'**  ^  .  J,  C*  -  S.  D 


Ajt  -  b3  !>a 


_  (u-x) 

^  }  “  +  }**■)  *  J* (i»i+  Jm)  Cj  -  Jv[)»i+~y*r)  Die. 


V-t'»4>0  &3  ♦  *>(>»*■ M  T-j>  (**+}*  3)  VjJ 

Boundary  conditions  at  the  bottom  (z  =  +  H^): 

_(*- y) 

Wso  o »■}■*■&)  A )  -  [t+Y'*)  &»*■  -V  C3  -  y.i>)  Jj«.  s 


If  J»>'  C3,=. 


If  *5  <1 


A3  *< 


*.  —  _______  {?>>«-  e  f  D^t 

2.  Resolution  of  the  system  of  10  boundary  equations. 

We  add  and  substract  the  surface  conditions 

2A,e  *  -1-  ^4  +  Ji)  C,eJ'*  -  [*-J') 

lb,  o  (^1-^0 

We  add  and  substract  the  first  two  conditions  at  the  first  interface 

2  A,  *  O f,)C,  *  *.e**  r  ]_**] 

z  b,i*  *■  {  A-i.)C,  *  (-»•*  ^  -  y 


A  4.8 


where 

[/U]  «  At.-f*-*1  5, 

We  replace  A^  and  Bj  by  their  values  in  function  of  and 

c,  £”y]  - »* 

0-r.)C,  »  (->*'•)  »,  *  [**>  ** 

We  solve  the  system 

c,.  {  [*»][».«. 

v  [2x,e*  -(>*x,)*x,x  v  Ias.)*1****']) J-^- 

n 

J),a  |  [A2][2Xie*  -  OX,)***  4 

+  ^2i\  «r’,ex,,f_  (A-hJ-,)  c  **  V  OX,)]  }  ^r- 

Vt«  -  O'O' l*2*  * 

We  transform  the  w-condition  utilizing  the  t  ^-conditions 

-  *,(•»,-.)  *,ZS'\  T,  [st(^0Cris,Cy-”\  D, 

and  with 

?t  «  F,  **  C"*^I 

S,  (w, —• ) 

C,.  3),^'%  ?a  [Czej'^-x;_  T>t] 

We  replace  and  by  their  values 

C*-  D,*-J‘,x= 

{l/OLC^O  -  0*K1X+  (><♦*) -i2**-  ^-xOe2-**^*] 

+  [ax^O-*7**)  *  2e2*e**_  2 ***]}. -±- 

»  1  !**]•*•  ♦  %- 
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Id 


[LA13  «■> 

4-  fit’ 

\  — 
j  V, 

LA,3  - 

nvi 

ft. 

lc* 

Rx 

R, 

Writing  ft** 

,  we 

obtain  the  system 

<-  *>a. 

4" 

£„-**)  Ct  + 

£>«+■  R^  )  = 

-  K 

V 

^  CiC*^’**  - 

3>i  =  o 

and  by  adding  and  substracting 

/H-rt»  *,)!>,  -  il 

Rl 

L^rt+R})  D,  -  A*. 

Hi 

We  add  and  substract  the  first  two  conditions  at  the  second  interface 
2^j  4-  ('U-Ji)Ci  *■  (h-Ji)  Sje  I/O 

2fcI«-£'<-’0 c.  +(-../■.)?, i^h-K 

where,  taken  in  account  the  bottom  conditions, 

lA%3-  L"+  4  ‘  J5i 

if 

I A»]  =  fc,  -r** ]  K 

We  replace  ^  ancl  ^  ^  their  values  from  the  first  interface  conditions 

C,[o*0  - 

v  J.  J_C -V»0 ]  =  y*5]  ♦  ^  Cl'»'x) 

C»  -  O J.-Rv)  e  M',")] 
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We  solve  the  system 

Ca»  |  [A>]£.2.Sie  e  ^  v  R»)  +  O  Jr+  »*] 

J),.  J  [Ai][2Sl 

v  -  {— joj>  v 

ni  1-  Vi 

7t.  U, (-r.-^)C2J,lv'"J 

for  m-  Oft  ,  'Vta  (>»-'*“*)  O'1  *  1*0 
We  transform  the  w-condition  utilizing  the  xrz-conditions 

C*.  -  rxC^t-”*)  J)i«-  ^  K  (*»>->}  *•  -  s»,[V) D  »^| 

and  with 

v, 

r*  (_m-0 

Cl  -  _  W*  T>»] 

We  replace  and  D2  ^  their  values 

Ct  -  Die^V"’0* 

-(x+j»  +  H)*c7'^'x)  -  O'-*-**)  *lty -*>~2*h-x>J 

=  {  (S,  ^  ®,  *j  -L-  =  k^c4;,,|s'v;-  t>»] 


with 
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$3= 

lA>3  -  Q»  t)>]  .  -  ih. 

This  relation  becomes 


If  **> 


A  +  e 


_x(z-y) 


]b,  ,  U+  hlmtiL ♦  4,1  P, 


)\  *><» 

It*  63 

+  \  -»  *•  X25^*'^  4.  (M  D3=  - 

L  J  <D, 

Together  with  the  r  ^-condition 

)[  r%>i 


HO 


]  ±y± 

*■ 


:,j  d 


we  solve  the  system 


H  H>» 


Ks.®i 


**£*'•»♦>*)  -lz“T) 

~7  ..;  ~T~  e  «•  -  H 


«k  «2^-. 


V3« 

C.'1 


A  4.12 


If 

v-| 


-2l*(7-y) 

*•  -J> 


jC*-v)  -j»U-v) 

Dv*  <&*  e  ”* 

6.  *T 

Va>  (^a,v  ±±i>  -f-V^-o 

(v>»*|^) 

3.  Relations  for  the  parameters. 

tx--  |  laOL2*«' t','x)-  (•,-., -no;’t''-V'(v-rt] 

A,»  *Jl»0 

The  relations  for  the  parameters  Ap  Bj,  Cj  and  Dj  are  given  in  §  1.3 
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Chapter  3.  The  four  layered  structure. 

1.  Supplementary  boundary  conditions. 

Boundary  conditions  at  the  third  interface  (z  =  +  Hg  +  H3): 

<r*:  A*  *  Cs  + 

Aa^K*J+  hA  *  Cz,  Jr*(Ka)  +  lA 

trx:  A,  - 

A*  •»-  J*  C*  e**  -  X-t, 

W  *.  (4*.^.*)  A*  -  I*  (*»%*■)»-*)  C*  -  *•►')  ^  * 

V.w  A*e^  ^  [A+-)**)  b*  V  fi*  (WavJaa)  Cax  ^  ”"  —  tt 

U_  :  Aj  V  C^+-V"^c^+  C"1*)^  V» 

Lw  Aac^”*\.  O*  ♦■)-*)  &a  v  (-4*^0  C* t,A  v-  iUj 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  +  H3  +  H^) : 

We  O 

+  Aa  Ja(»u*>a)Ca  -**  (►u^HOJUfc  ^  e« 


'{ 

*A>  » 

Ca  "  o 

A<iW. 

«.  *1  ^ jlV-e) 

'1 

5“  A  <) 

A  a  3  0 

va  v 

.  w> .  fc4e-tM’:'uv•*,  - 

*■>**) 

A4.  14 


2.  Expression  of  the  boundary  conditions  in  matrixform. 

2.1.  At  the  third  Interface. 

In  the  equations  at  the  third  interface,  A^,  or  C^,  is  replaced  by  its 
value  obtained  from  the  fixed  bottom  condition. 

We  write  the  conditions  at  the  third  interface  in  matrixform 


»<0T 


We  invert  M, 


s3  -Ja  \ 

A  ^  **(”**)*- ^  -s»  J*  ^  -S>>S2~y}  ] 

-*>(*>+ H)  (-**►*)  s>  . 

»»M  ,  .HM  M**v)  I 

s4>  1 


-70-*  ,  +  £4iv^:L  n 

1*V«  ,, . 

umi™**]  u  c^>  t- 


<‘ 

AWiA+)ia) 


-2*4 1  Kft) 
ne 

-2X4(HjO 

-;<k 


^vjA)lKotn  L*|vu,*H[->,  +  e'  ^  J 

I  (>ivUa1  \-  1  *■  -* 

UA  *  '  1  f  -^*,)(S*>n  .  ,  xr  -ZUMI 
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We  write  A^,  B3,  C3  and  D3  In  function  of  and  D^, 


A  ♦  -2?  0*)P* 

/»3  *  — - 


*3  = 


2  tty  ftp.  Ht £h).S>a  »  ±  (a'y) 


2^-*  a) 

C  McOQ.  €m  v  ^  NgfrO-  fl<0)D4 

*  2*>  £*-*>) 

*  ih»(40-  Wc[u).b4  * 

l)  5  x  " 

where  Mg(1,j)  are  the  constants  in  M^. 

We  write 

TjH  -  X-H'ClVMCO*) 

?j»  «  .£  Msfi’O  M 1 62) 

so  that 


T« 

^.2  . 

T„^ 

Pnet2'y) 

?3. 

V 

•0 

W«- 

<0  *»(*•*) 

"43  e 
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2.2.  At  the  surface 

Adding  and  substracting  the  surface  conditions  we  obtain 

7  b,  .  >l_  (—■ r.)  C. .-M- (**'•)*' 

2.3.  At  the  first  interface. 

We  add  and  substract  the  first  two  conditions 

lA.  *  (~»)C,  ,  ■  *’cty"‘K  ^ 

■  LA*] 

2b ,C%  ♦  H*-J0  lA23 

We  replace  Aj  and  Bj  by  their  values  obtained  from  the  surface  conditions 

+  (■*-(.)!, I*'"'-  e.“J  =  I; A 73  -  t* 

=  l**]-*' 

We  solve  the  system 

allA2]l2r,^^,c‘  (*+*)£*<■ 

l  _'!*  -x.in  \  4 

¥  YJ25-, cT7*  -  4  {*-*')*-  e  -J  )•  *^5, 

The  positive  exponent  ex  has  disappeared. 

3),.  h*0  .  (*■>'-'!  «**  *  ^ 

*■  [u.,-’iw  -  V  ~vT 

For  V>1  =  «3  7>)  =• 


M.  17 


We  transform  the  w-cogdition  utilizing  the  x  ^-conditions 
r,C»n-i)  c,  -  C*'~ 0  D>e^'*=  T, 

Ft 

Tt), 

We  replace  and  Dj  by  their  values 
C,-  D,«Er,x= 

For  wi=  oC>  H-i  =  (>'-S0  >  ° 


*R’V  vT= 


vcte  '-'r’2-  fp3 


Writing  R^»  ^3l  ,  we  obtain  the  system 

f?* 

6»  *  C.^’v  J,  -  *>T>^  -■£ 

4 ucjiti-*- 

and  by  adding  and  substracting 

lA,ct7-^=.  *»)•>»- 


* 
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2.4.  At  the  second  interface. 

We  add  and  substract  the  first  two  conditions 

2Aa  4-  [*+Jt )C%  *-  Vi*311'*"* I  A^'yl  fc*  +■  Cvlf*(a'y)+ 

-  1**1 

abl£^"XV  ^-rx)Ca  V  [A  3] 

We  replace  A2  and  B2  by  their  values  obtained  from  the  first  interface 
conditions 

♦  CA33  +  I: 

CZ]_^U)  e  J 

v  D2[c^+x0^ty’X-  +  X)]  -  ^A^  +  TT^(y 

We  solve  the  system 

s.  L'  yiH>.  v2 


«,  U  C-»*  ^ 

The  positive  exponent  e^y-x^  has  disappeared. 

Da-  ^  l**] Lax,^'f"y)_  i****- *0**^"*^  ,y  *■  1  J 

Tor  vm*o»  ,  Vt- 
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We  transform  the  w-condition  utilizing  the  1  ^-conditions 

W  =  k,  int^r 

S-i  Oa“') 

Cx_  XcC^”*-*** 

We  replace  C2  and  D2  by  their  values 

CX-W,C*"5* 

1  Lai3[C^->  ri-sO«-  ^ 

’  ^  -2N-X)  _2J-,|y-x')-i 

W"-  ai^V^  -^r 

s.L 

«  ^[A^Q,  -  ®A  ^ 

Tor  „-«  ,  fi>,-  C>l-*»*SO  , 

Writing's  \£^  ,  we  obtain  the  system 

A b3  V  -®>C^’>|,iy)<.(!I!,1>j=- ^ 

A  ®3  *  j.CiC’W*^  .  I»t>»  -« 


and  by  adding  and  substractinq 

4  (,+  rs4>)  CSC,,1’'"V)  »  t-  s*  *-  '  $7 
11b,  4  *•  '  77 
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3.  Resolution  of  the  system  of  boundary  conditions. 

We  have  from  the  boundary  conditions  at  the  third  interface  (  §  2.1) 

A^a  - — -  V  Pj|.  ^>4  + 

d  -  A  [  ^a3>*] 

3 

X)3.  -  +  Pax  tuj  eJ>^"V) 

We  replace  /*3,  By  C3  and  D3  in  the  last  equations  of  §  2.4 

eT^*'v^-*-  £a*  *)-&■>)  ?>«  V^+  (*-5*+  ®i)  PA»e-^  Ba 

»  (W>  -<s>)  f«c,,|’'y>»  (*>.**,)I!«1,>|m1] 

l-LV^-'h  .  ^f**»0&*,*M}»* 


-  -7^(A-Yi)  ^iL 


We  solve  the  system  in  and 


i^\^V3 
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r 


»r» ^ 

6,  L 

v  >,  ,,*„£(*-»)-].  -I- 

J  Vs 


The  positive  exponents  e^2"^,  es3^2"^  and  e^2-^ *es3^z'-y^  have 
disappeared. 


D*  =  2**  (*-*>)  4i-  I  Pa,  £M“y)-  **, 
to.  L 


V^-  d]_CP«  P^x-^p20  +  y)e Y) 

¥  (»4  ^2fz'y‘> 

,  _  s  -lsJz~y) 

*  ^  +  /*-  l&3)  (?i<^  -  HzP^O  e 

+  (^-h  .ii)[^h-,!‘^ 

V  ™  =  *  7,  «  1>-M  ♦  ( h* ?2"  -  721  Rt2) 

The  numerators  of  8^  and  tend  both  to  zero  and  the  denominator 
tends  to  a  constant  value. 
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4.  Values  of  the  parameters  A^ ,  . 

4.1.  Values  of  the  parameters  A3,  B3,  C3,  D3> 

The  values  of  the  parameters  A3,  B3,  C3.and  D3  are  obtained  from  the 
boundary  conditions  at  the  third  interface  in  which  B^  and 
are  replaced  by  their  values  from  §  3. 

+  p,aD^  ■  rVr 7 
a  U  1*1  (/*-**) 

ll*  fpz,  fix 

L  2 Si  [*-**) 


-  f  (P,!.?!,  -  P»  P^)  '  «- 

7 o  .p»i -  fcv *•) i ”4|,T'  *  M*- -  *■  *0* * y,J  ^ 

.  P21-  ^ 

The  positive  exponent  e^z_y^  has  disappeared.  Although  the  presence 

of  the  constant  S3(P22^41  "  P21^42^  t^ie  numerator  converges  to  zero 
because  of  the  factor  Qj. 


Cj*  j^Pv  TSi  ^  A J 


2  i>  (*->'*) 


%.!  cv*>  -  iv».  **‘M*  if ».v 

CD,  L 

4-  -  Pva*')  «7Z  VM  — - 

v  -1  T7, 


t>3*  [_  ?AY  *  ^1.0*3 

2/j 

«  T  I  Va\  .  Pa  -  V  ?4*3  *  2  **  *  (?*»•  h*  -  ?4  »•  p“) 

L 

4.2.  Values  of  the  parameters  A2,  B2,  C2,  D2. 

The  values  of  C2  and  D2  are  obtained  from  the  relations  established  in  §  2.4. 

Ci.  i— !  -  (-*.*  *oj. 

1  1  \A»tl*‘v)  *  B,  *  CiC'^-’h  t>»] 

4-  ^-.vu)  £’,ll'-X)-  1*«  I'H  -  (-J.)  ^ 

Dx=  *—  j 

*■  S>5  ►c>c,»l1'y)  ♦  i>3] 

.■  h.  (»-*>)]  ] 

The  numerator  in  D2  converges  because  of  the  presence  of  the  factor  R2> 

The  value  of  A2  is  obtained  from  the  relation  established  in  §  2.4. 

Ax,  i.  ^ A,  C (’■-*>  ,  s,  *  C1<r,*lt->'\'D3 

The  value  of  B2  is  obtained  from  the  last  relation  of  §  2.3. 

-u-*>')ClCilL't-*K  (44.x^R0T)2  ■* 
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4.3.  Values  of  the  parameters  A^,  D^. 

The  values  of  the  parameters  A^  and  are  obtained  from  the  relations 
established  in  §  2.3. 

C,  ■=  — —  |  ^  ■ 

[A.ilv-J  .  3>,] 

M  3  y 

For  the  determination  of  A^,  we  need  the  value  of  D^e  1  . 

-b,  £r'x=  A__  j  ^  c'C*’’  -  <■  [a~>0 

[A itb-*)  v  a,  <-  Cjt’ilv-1'1  v 

*-  [lr,Cx  e™*. 

A,= B2  V  Ilj  -  [*■>?, )C, 

The  values  of  and  are  obtained  from  the  surface  conditions. 

B,=  — — ■  W,  -  (.A.foA.e'  -2f,c, 

"J)  -  -  *  -[aiv5'i)Ci£  1 

1  ~  s,-' l  L 

5.  Relations  for  the  stresses  and  the  displacements. 

The  relations  for  the  stresses  and  the  displacements  are  completely 
the  same  as  those  developped  in  appendix  3,  by  replacing  the  parameters 
A^ ,  D.  by  their  adequate  values. 

The  relation  for  the  verticale  displacement  is  again  undeterminated 
at  the  origin  (  m  =  0).  The  problem  is  solved  in  exactly  the  same 
way  as  developped  in  appendix  2  (§  4.). 
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* 
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*************************************************************************** 

Four-layered  System  Program  : 


This  program  is  available  in  two  versions: 

-  EXECUTABLE  version,  in  which  the  executable  program 
is  made  up  of  only  one  block  with  automatic  loading. 

-  SOURCE  version,  in  which  all  controls,  data,  and 
all  modules  are  in  separated  tiles. 

Floppy  disk  contents  : 


In  -following  text,  x  means  A  =  ANISOTROPIC 

I  =  ISOTROPIC 
y  means  F'  =  PARTIAL 
S  =  SLIP 

Each  -floppy  disk  contains  following  main  files: 

AUTOEXEC.BAT  - >  automatic  program  loading 

LOGQ.E-tAT  - >  display  of  introduction  logo 

FLxyLQ.TXT  - >  introduction  text  on  screen 

FLxyNO.LOG  - >  introduction  text  for  printer 

FLxyNO.TXT  - >  this  notice 

♦EXECUTABLE  floppy  disk 

This  disk  contains  the  following  files  in  addition  to 
the  main  files: 

FLxy.EXE  - >  executable  program 

FLxy.DAT  - >  data  file  for  demonstration 

FLxy.LST  - >  result  file  for  demonstration 

♦SOURCE  floppy  disk 

This  floppy  disk  contains  the  following  files  in  addition  to 
the  main  files: 

FLxy.VER  - >  revision 

FLxyl.FOR  - >  main  module 

FLxyZ.FuR  - >  subroutines  DOMEC 

DOTRA 

FOCAL 

POINT 

PAS 

CHECH 

ERROR 

FLxy3.F0R  - >  subroutines  ECHDE 

VINIT 

ZERO 

FINIT 

FLxy4.F0R  - >  subroutines  F'4442 

PCT22 

PCT42 

S0M42 

S0M22 

F’4444 

P2442 

CONST 
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FLxy5. FOR 

- >  subroutines 

BESJ1 

BJ0J2 

SURFA 

C0UC1 

C0UC2 

C0UC3 

CQUC4 

FLxyS. FOR 

- >  subroutines 

TITRE 
MOD  IF 
SIMPS 

FLxy7. FOR 

- >  subroutines 

COMSU 

FONC 

CC0U1 

FONC  1 

F0NC2 

F0NC3 

F0NC4 

F0NC5 

F0NC6 

ECHEF 

FLxyS. FOR 

- >  subroutines 

IMDON 

IMRES 

FLxy9. FOR 

- >  subroutines 

MENU 

ECDON 

LEDON 

AFDMC 

AFDTR 

AFPOC 

Configuration: 


IBM-PC  (G,XT,AT)  with  at  least  256KB,  1  ou  2  diskette  drives, 

80  col.  screen  (monochrome  or  color),  math,  coprocessor, 
matrix  printer. 

Running  of  the  program: 


In  the  next  presentation  < ENTER >  means  action  of  key  < 
Insert  the  EXECUTABLE  disk  into  drive  B. 

Load  DOS,  if  necessary,  then  type  DIR  B:  and  PATH  A:  \ 
FLB>  FLxy  < ENTER > 

Printing  of  the  results: 


on  screen  :  FLB>  type  FLxy.LST  <ENTER> 

on  printer  :  FLB>  type  FLxy.LST  >prn  <ENTER> 

Consultation  of  this  notice: 


on  screen 
on  printer 


FLB>  type  FLxyN0.TXT  < ENTER > 

FLB>  type  FLxyN0.TXT  >prn  < ENTER > 
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*************************************************************************** 

Disk  preparation  -for  normal  running: 

Format  a  system  disk  with  CQMMAND.COM,  and  tiles  needed  -for 
running  AUTOEXEC.BAT. 

Preparation  of  CONFIG.SYS  with  files=10,  device=ansi,  buffers=10. 
Alterations  in  source  files: 


By  text  editor  EDLIN  or  any  other  text  editor,  program 
statements,  may  be  altered. 

Commands  for  EDLIN  are: 

FLB>  edlin  xxxxn.for  <ENTER> 

where  - >  file  name  to  be  modified  (xxxxn.for) 

nD  — >  erase  line  number  n 

n  — >  displays  line  n  for  alteration,  type  the 

correct  statement. 

nl,n2L  — >  displays  lines  between  number  nl  and  number  n2 
E  — >  ends  session  and  returns  to  DOS 

Compi 1 ati on : 


After  alteration,  the  new  version  of  the  module  has  to  be 
compi 1 ed . 

Insert  Professional  fortran  compiler  into  drive  A. 

FLA>  b:  < ENTER > 

FLB>  path  a: \  < ENTER > 

FLB>  profort  xxxx/li  < ENTER > 

Linkage: 


After  correct  compilation  has  taken  place,  a  new  executable 
program  has  to  be  created. 

Insert  diskette  with  FORTRAN  libraries  into  drive  A. 

FLB>  link  FLxy l+FLxy2+. FLxy, CON: ;  < ENTER > 

Note  : 


To  obtain  introduction  logo  on  printer 
FLB>  type  FLxyNO.LOG  >PRN  < ENTER > 
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*************************************************************************** 

Presentation  of  different  possible  screens: 


FD-USAE-vers.  2.00  1986 

MAIN  MENU  scren  : 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 

******************  a######*####*##*####*#*#*##**##*###*##*****##*#*##**  #:*#***** 


1 . 


4. 

5 . 

6 . 
7. 


Data  retrieval  in  a  Pile 
Data  saving  in  a  file 


Screen 

displaying 

and/or 

al ter ati 

on 

o-f 

Screen 

displaying 

and/or 

al terati 

on 

of 

Screen 

di spl ayi ng 

and/or 

al ter ati 

on 

of 

Input  of  intermediate  depths 
Program  start 


system  data 

traffic  data 

computati on  coordi nates 


Your  choice 


0  for  stop  —  : 


Screen  of  CHOICE  1: 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 

****************************************************************************** 
Data  retrieval  in  a  file. 


A  FLxy.DAT  file  contains  base  data. 

The  user  can  define  another  file  whose  name  nas  to  be  written  in 
•8  characters  (format  XXXX.DAT). 

Name  of  chosen  file  or  FLxy.DAT  : 
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♦  **##****. **♦***!!<*■  **t:*.****)*:*)|<******^^*  ******************  Y^*|!!(<***))<**Y**^***^* 

Screen  oT  CHOICE  2: 


Strains  and  stresses 

in  USAE— 2 

a  -four-layered  system. 

^  ^  ^  ^  ^  4'  ^  ^  ^  ^  ^  ^  ^  ^  ^  X  X  dr  X  d  \L  X'  d>  X  X  X  rV  X  X  X  X  X-  X  X  X  X  x  x  x  x  x  x,  X  x  x  x  x  X  x 

*  *  *  *  ^  t*****^*^1  t  ^  ^  *  m  *  m  *  t  *  m  *  t*  *  m  m  *  *  ■*  *  m  *  *  *  m  m  *  *  *  ^  *  *  m  *  ^  ^  ^  ^  t  ^  ^  ^  f  ^  ^  ^  ^  ^  m  t  ^  ^  m  t- 

Data  saving  in  a  tile. 


A  FLxy.DAT  tile  contains  base  data. 

The  user  can  detine  another  tile  whose  name  has  to  be  written  in 
3  characters  < -format  XXXX.DAT). 

Name  ot  chosen  tile  or  FLxy.DAT  : 


.  Screen  ot 


CHOICE  3: 


Strains  and  stresses 

in  USAE -2 

a  tour-layered  system. 

Jit#####*###########*:*##*###**####*#*#*######*###  ##*#*#*#.  %  if:#*#*  #:****:**  *** 


System  data. 


Layer 

Modul us 

Poi sson ’ s  r . 

Thickness 

1 . 

400000. 0 

0.  16 

20 . 0 

9 

1 UOOOO . 0 

0.25 

20.0 

3  • 

10000.0 

0 .  50 

30.  0 

4. 

1000.0 

0.50 

9000. 0 

Friction  ra. 


1 . 0000 
1 . 0000 
1 . 0000 


•For  return  =1 


or 


alter 


=2 


ISIC  FOUR-LAYERED  (FL:;y>  USAE  Page  6 

*************************************************************************** 


Screen  ot  CHOICE  4: 


Strains  and  stresses 

in  USAE-2 

a  tour-layered  system. 

^  d>  d-  d  4  ^  d  d  d  d  d  d  d  d  d  ddddddddddddd  d  dddddddd1 

^  f  ^  f  f  ^  f  ^  ^  ^  ^  m  d  m  ^  ®  ^  ^  m  ^  m  m  m  ^  ^  ^  m  m  m  ^  ^  m  *  ^  ^  m  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *•  ^  t  / 

Number  ot  circular  loads  =  2 

Radius  Pressure  Dist.  ;•:  Dist.  y 


1 1 . 450  7 . 900  0 . 000  0 . 000 

1 1 . 450  7 . 900  34 . 350  0 . 000 


For  continuation  =1  or  alter.  =  2 


Screen  at  CHOICE  5: 


BAE  —  A 

*  *  *  *  *  ******** 

Number  ot  computation  coordinates  =  2 

0.000  0.000 

17.175  0.000 

For  conti nuati on=  1  or  alter.  =  2  -  : 


Strains  and  stresses 

in  Ui 

a  tour-layered  system. 

*****************************************************************: 


Screen  ot  CHOICE  6: 


Strains  and  stresses 

in  USAE-2 

a  tour-layered  system. 

****************************************************************************** 

Positions  ot  stress  computations  in  depth  out  ot  intertaces. 

(max  22) 
number  ot  positions  : 
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a************************************************************************** 

Screen  of  CHOICE  7: 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 

X  ^  ^  ^  ^  .1.  ^  ^  ^  ^  ^  ^  yL  X  tL  ^  ^  ,1  ^  ^  ^  - 1 .  ^  ^  ill  ^  ^  >lL  ■  i »  )!•  tV  dr  ^  d#  df  d1  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d  d 

t  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  m  m  *  *  *  *  T  *  *  ^  *  *  T  *  *  ^  ^  ^  *  *  T  *  *  ^  ^  ^  ^  -T  t*  *  *  <f  T  *  ^  ^  ^  ^  ^  ^  ^  ^  • 


initial  computation  interval  <0.1  is  generaly  small  enough) 


Next  screen: 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 


choice  of  scale  : 

allowed  choices  : 


1  - 

thickness 

o-f 

lrst  layer 

'“>  _ _ 

thickness 

of 

2  first  layers 

3  - 

thi ckness 

of 

3  first  layers 

4  - 

thickness 

of 

4  layers 

5 

load  radi1 

LIS 

suggested 

soluti on 

:  5 

your  choice  : 


Screen  for  execution: 


Strains  and  stresses 
i  n 


USAE-2 


a  four-layered  system. 

’***M:M**********************:M*****iM**************************************** 


computation  start .  be  patient!! 

m  =  0. 10 

be  even  more  patient. . .  !  ! 


1 
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**************************>:************************************************ 

Screen  for  results: 


Strains  and  stresses 

in  USAE-2 

a  four— layered  system. 

****************************************************************************** 


A  file  FLxy.LST  contains  base  results. 

The  user  can  define  another  file  whose  name  has  to  be  written  in 
8  characters  (format  XXXX.DAT). 


Name  of  chosen  file  or  FLxy.LST 
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IK*****#*********#**********#*****************************)*!*#**.  ************* 

Sample  o-f  results  (FL;:y.LST) 


Four-layered  system  program  isotropic  -  partial  -friction 

USAE-2  departeient  des  constructions 

i  sic 

av.de  l’hopital,  27  h 
7000  ions  belgi que 

mmmmtitmutttmtttmmtiiitmtmmummmtmttttmttmiuttmmmtttuummuttttttiu 

Mechanical  data 
*************** 

Young's  modulus  P.’s  ratio  thickness  friction  r 


4000 00 . o 

0.  16 

20 . 0u0 

i . 

0f» 

1 00000 . 0 

0 . 25 

20. 000 

i . 

00 

1 0000 . 0 

0.  50 

30 . 000 

i . 

0'J 

1 000. 0 

0 .  50 

9000 . 000 

traffic  data 

************ 

v 

0.000  0.000 

1 1 . 450  7 . 900  34 . 350  0 . 000 


load  radius  pressure 

1  11.450  7.900 
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tiittttttttttttttttttitttttttttttttttttttttttitttitttuttttttttitimttim 


position 

1  X 

=  0.000 

y 1 

0.000 

depth 

layer 

sx 

S2 

ty: 

txz 

txy 

0.000 

1  tt 

10.1600 

11.6793 

7.9000 

0.0000 

0.0000 

0.0000 

20.000 

1  tt 

-3.4792 

-4.5018 

1.4988 

0.0000 

-0.4167 

0.0000 

20.000 

2  tt 

-0.5790 

-0.3163 

1.4989 

0.0000 

-0.4166 

0.0000 

40.000 

2  tt 

-2.2733 

-2.5873 

0.1873 

0.0000 

-0.0611 

0.0000 

40.000 

3  tt 

-0.173a 

-0.1997 

0.1873 

0.0000 

-0.0611 

0.0000 

70.000 

3  tt 

-0.3907 

-0.4035 

0.0540 

0.0000 

-0.0046 

0. 0000 

70.000 

4  tt 

0.0095 

0.0082 

0.0540 

0.0000 

-0.0046 

0. 0000 

9070.000 

4  tt 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

depth 

layer 

Si 

s2 

s3 

epsl 

eps2 

eps3 

0.000 

1  tt 

11.6793 

10.1660 

7.9000 

U.2197E-04 

0.1758E-04 

0.1 101E-04 

20.000 

1  tt 

1.5335 

-3.5138 

-4.5018 

0.7040E-05 

-.7597E-05 

-.1046E-04 

20.000 

2  tt 

1.5793 

-0.6595 

-0.8163 

0.194BE-04 

-. 8502E-05 

-.1Q46E-04 

40.000 

2  tt 

0.1838 

-2.2749 

-2.5873 

0.1 404E-04 

-.1675E-04 

-.2066E-04 

40.000 

3  tt 

0.1973 

-0.1836 

-0.1997 

0.3890E-04 

1924E-04 

-.2066E-04 

70.000 

3  tt 

0.0540 

-0.3908 

-0.4035 

0.4512E-04 

-.2160E-04 

-.2352E-04 

70.000 

4  tt 

0.0545 

0.0091 

0.0082 

0.4582E-04 

-.2230E-04 

-.2352E-04 

9070.000 

4  tt 

0.0000 

0.0000 

0.0000 

0. OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

depth  1 

layer  u(x) 

v  (y ) 

•t  <r ) 

ex 

ey 

e: 

0.000 

1  tt  -.2184E-03 

-. 1909E-10 

-.9046E-02 

0.1758E-04 

0.2197E-04 

0. 1101E-04 

20.000 

1  tl  0.9985E-04 

0.8642E-11 

-.8835E-02 

-.7497E-05 

-. 1046E-04 

0.6939E-05 

20.000 

2  II  0.9885E-04 

0.8642E-11 

-.8835E-02 

-.7497E-05 

-.1046E-04 

0. 1848E-04 

40.000 

2  It  0.3017E-03 

0.2637E-10 

-.8561E-02 

-. 1673E-04 

-.2066E-04 

0.1402E-04 

40.000 

3  II  0.3017E-03 

0.2637E-10 

-.8561E-02 

-. 1673E-04 

-.2066E-04 

0.3739E-04 

70.000 

3  It  0.3844E-03 

0.3360E-10 

-.7432E-02 

-.2159E-04 

-.2352E-04 

0.4511E-04 

70.000 

4  tl  0.3844E-03 

0.3360E-10 

-.7432E-02 

-.2159E-04 

-.2352E-04 

0.4511E-04 

9070.000 

4  tt  O.OOOOE+OO 

O.OQOQE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 
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timttttmtittiimtiittttmtttstimttttitttttttitttttmttitmttttttt 


position 

2  x  ' 

=  17.175 

y  = 

0.000 

depth 

layer 

sx 

sy 

S2 

tyz 

tX2 

txy 

0.000 

1  II 

3.0499 

8.8259 

0.0000 

0.0000 

0.0000 

0.0000 

20.000 

1  II 

-2.0893 

-4.3575 

1.1255 

0.0000 

0.0000 

0.0000 

20.000 

2  tl 

-0.3243 

-0.8505 

1.1257 

0.0000 

0.0000 

0.0000 

40.000 

2  II 

-2.4172 

-2.7155 

0.1997 

0.0000 

0.0000 

0.0000 

40.000 

3  II 

-0.1828 

-0.2077 

0.1997 

0.0000 

0.0000 

0.0000 

70.000 

3  It 

-0.4084 

-0.4166 

0.0554 

0.0000 

0.0000 

0.0000 

70.000 

4  II 

0.0090 

0.0082 

0.0554 

0.0000 

0.0000 

0.0000 

9070.000 

4  II 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

depth 

layer 

si 

s2 

s3 

epsl 

eps2 

eps3 

0.000 

1  II 

8.8259 

3.0499 

0.0000 

0.2084E-04 

0.4094E-05 

-.4750E-05 

20.000 

1  tl 

1.1255 

-2.0893 

-4.3575 

0.5392E-05 

-.3930E-05 

-.  1051E-04 

20.000 

2  II 

1.1257 

-0.3243 

-0.8505 

0. 1419E-04 

-.3931E-05 

-. 1051E-04 

40.000 

2  II 

0.1997 

-2.4172 

-2.7155 

0. 1483E-04 

-.  178BE-04 

-.2161E-04 

40.000 

3  II 

0.1997 

-0.1828 

-0.2077 

0.3949E-04 

-. 1788E-04 

-.2161E-04 

70.000 

3  II 

0.0554 

-0.4084 

-0.4166 

0.4679E-04 

-.2278E-04 

-.24O1E-04 

70.000 

4  II 

0.0554 

0.0090 

0.0082 

0.4679E-04 

-.  227BE-04 

-.2401E-04 

9070.000 

4  it 

0.0000 

0.0000 

0.0000 

O.OOOOE+OO 

0.0000E+00 

O.OOOOE+OO 

depth  layer 

utx) 

v  (y  J 

*(:) 

ex 

ey 

tl 

0.000 

1 

II 

0.2272E-18 

-. 1565E-10 

-.8979E-02 

0.4094E-05 

0.2084E-04 

-. 4750E-05 

20.000 

1 

II 

-. 4354E-19 

0.7889E-11 

-.8961E-02 

-.3930E-05 

-.  1051E-04 

0.5392E-05 

20.000 

2 

II 

-.3448E-18 

0.7889E-11 

-.8961E-02 

-.3931E-05 

-. 1051E-04 

0. 1419E-04 

40.000 

2 

It 

-.7092E-18 

0.1622E-10 

-.8698E-0? 

-.1788E-04 

-.2161E-04 

0.1483E-04 

40.000 

y 

J 

II 

-.7092E-18 

0. 1622E-10 

-.8698E-02 

-.178BE-04 

-. 2161E-04 

0.3949E-04 

70.000 

3 

It 

-.7879E-18 

0.1803E-10 

-.7520E-02 

-.2278E-04 

-.  2401E-04 

0.4679E-04 

70.000 

4 

II 

-.7879E-18 

0. 1803E-10 

-.7520E-02 

-.  2278E-04 

-.2401E-04 

0.4679E-04 

9070.000 

4 

It 

0 .  OOOOEe-OO 

O.OOOOE+OO 

0.00008*00 

O.OOOOEtOO 

O.OOOOEtOO 

0. 0000E+00 

i 
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Utilized  symbols: 


sx 

sy 

sz 

tyz 

txz 

txy 

si 

s2 

s3 

epsl 

eps2 

eps3 

u(x) 

v(y) 

w(z) 

ex 

ey 

ez 


normal  stress  in  the  x-di recti  on 

normal  stress  in  the  y-direction 

normal  stress  in  the  z-direction 

shear  stress  in  the  yz  plane,  parallel  to  y  or  z 

shear  stress  in  the  xz  plane,  parallel  to  x  or  z 

shear  stress  in  the  xy  plane,  parallel  to  x  or  y 

maximum  principal  stress 

medium  principal  stress 

minimum  principal  stress 

principal  strain 

principal  strain 

principal  strain 

displacement  in  the  x-di recti  on 

displacement  in  the  y-direction 

displacement  in  the  z-direction 

strain  in  the  x-di recti  on 

strain  in  the  y-direction 

strain  in  the  z-direction 


The  normal  stresses  are  taken  positive  when  they  produce  compression 
and  negative  when  they  produce  tension. 


